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"What makes you think you can discover anything? Who are you?" 
"Nobody. Nobody at all. But the secrets of the Universe don't mind. 
They reveal themselves to nobodies who care." 

- The Outer Limits: The Galaxy Being 
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Summary 

In this thesis, we use the motion of the Local Group of galaxies (LG) through the 
Universe to measure the cosmological parameter of non-relativistic matter density, Q^. 
For that purpose, we compare the peculiar velocity of the LG with its gravitational 
acceleration. The former is known from the dipole of the cosmic microwave background 
radiation and the latter is estimated here from the clustering dipole of galaxies in 
the Two Micron All Sky Survey (2MASS) Extended Source Catalog. We start by 
presenting the general framework of perturbation theory of gravitational instability in 
the expanding Universe and how it applies to the peculiar motion of the LG. Next, we 
study a particular effect for the dipole measurement, related to the fact that a nearby 
Local Void is partially hidden behind our Galaxy. We then describe in detail how 
we handled the 2MASS extragalactic data for the purpose of our analysis. Finally, we 
present two methods to estimate the density Q^a, combined with the linear biasing b into 
the parameter /3 = Q'^^ /b, from the comparison of the LG velocity and acceleration. 
The first approach is to study the growth of the 2MASS clustering dipole with increased 
depth of the sample and compare it with theoretical expectations. The second is 
to apply the maximum-likelihood method in order to improve the precision of the 
measurement. With both these methods we find /3 ~ 0.4 and ^ 0.2, which is 
consistent with various independent estimates. We also briefly mention some future 
prospects in the field. 



This Thesis presents and expands the results included in the following publications: 
BiLiCKi & Chodorowski (2010), MNRAS, 406, 1358 - Chapter 2, 
BiLiCKi ET AL. (2011), ApJ, 741, 31 - Chapters 3 and 4. 

Additionally, in Chapters 4 o-nd 5 we apply the findings of Chodorowski et al. (2008), 
MNRAS, 389, 717. 

Material unpublished at the time of Thesis submission is contained in Chapter 5. 
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1 

Introduction 



bservational cosmology has been developing rapidly in the recent years. Progress 
in detector technology, usage of still bigger telescopes in more and more astro- 
nomically favorable places (including Antarctica and space), supported by numerous 
theoretical investigations and computer simulations, resulted in the emergence of a 
model of the Universe, which is considered to be standard: the Lambda-Cold Dark 
Matter cosmological model (ACDM). It is based on the cosmological principle: the 
Universe is homogeneous and isotropic when smoothed over large scales. This assump- 
tion is supported by various observations, including these of the cosmic microwave 
background (CMB) and the large-scale distribution of galaxies. 

Among the great successes of the ACDM model, one should mention its ability 
to perfectly reproduce the observed properties of the CMB. The angular power spec- 
trum of the temperature distribution of this radiation is now accurately measured up 
to the multipole moment of £ ^ 3000, thanks to several years of observations of the 
WMAP satellite (Jarosik et al., 2011), as well as new ground-based facilities, such 
as the Atacama Cosmology Telescope (Das et al., 2011) or the South Pole Tele- 
scope (Keisler et al., 2011). These measurements, as well as those of the CMB 
temperature-polarization power spectrum (currently known up to £ ^ 700, Jarosik 
ET AL. 2011), are perfectly fitted by a simple spatially flat model with non-zero cosmo- 
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logical constant A. Another evidence strongly supporting the validity of this model on 
the largest scales are numerous measurements of baryon acoustic oscillations (BAO), 
i.e. patterns of baryonic matter clustering at certain length scales due to acoustic waves 
which propagated in the early universe. BAO give us a 'standard ruler', whose size 
measured at various cosmological epochs is consistent with predictions of ACDM (e.g. 
Percival et al. 2010; Blake et al. 2011; Beutler et al. 2011). Last but 
not least, observations of 'standard candles' such as supernovae la, which actually gave 
first firm arguments for a non-zero value of A (RiESS ET AL., 1998; Perlmutter 
ET AL., 1999), additionally confirm that this model is a very successful general de- 
scription of our Universe (e.g. Kessler et al. 2009). 

On smaller scales we cannot of course treat the Universe as perfectly homogeneous, 
because galaxies not only do exist, but also are not arranged randomly and form struc- 
tures such as clusters, superclusters, filaments or walls, separated by huge voids. For 
that reason, there are claims that the observed inhomogeneities cannot be neglected 
when analyzing the big picture of the Cosmos, for instance by accounting for the 
so-called backreaction (see e.g. Buchert 2011 and references therein). However, in 
general these effects are believed to be weak and owing to the considerations presented 
in the previous paragraph we will not address them here. On the contrary, from now 
on we will assume that the cosmological background is statistically homogeneous and 
isotropic over large scales. What is more, it is known that gravity is the major force 
responsible for large-scale structure evolution - the other three fundamental interac- 
tions are negligible on cosmic scales. This means altogether that the Universe as a 
whole can be described in the framework of Friedman- Lemaitre models, obtained by 
solving the Einstein equations under these assumptions, and any perturbations with 
respect to this idealized picture will be considered within this model as a background. 



1.1. FRIEDMAN-LEMAITRE MODELS 



1.1 Friedman- LemaTtre models 



The cosmological principle expressed mathematically says that the metric of the 
Universe (the 'background') is given by the Robertson- Walker form: 



+ x' dr + sin"^d0^ 



1 — kx"^ 

(Robertson, 1929; Walker, 1933), where a{t) is the scale factor, t is the cosmic 
time, k is the curvature constant and {x, 6, (p) are spatial (spherical polar) coordinates. 

Another assumption underlying the relativistic cosmology is that the energy-momentum 
tensor of the matter field, T^,^, takes the form proper for the ideal fluid: 

T^u = (p + p)u^u^ + pg^u , 

where is the four-velocity of matter, p stands for energy density, p is pressure and 
g^i, the metric tensor. 

Writing down the general form of the Einstein equations (Einstein, 1916), 

(with R^y the Ricci tensor and R the Ricci scalar) and plugging Eqs. (1.1) and (1.2) 
into it, we get the following two independent Friedman- Lemaitre equations: 

tt2 f(^V k A ^— ^ 

^ = -^(p + 3p)4 (ID 

a 6 6 

(Friedman, 1922; Lemaitre, 1927), where H {Hubble parameter), a, p and p are 
functions of the cosmic time. 

The major components of the energy-density budget of the Universe can be treated 
as perfect fluids with the equation of state (EOS), p{p), given by a simple relation 



p = wp . 



[16] 



In particular, for the purpose of cosmological analyses the pressure of non-relativistic 
matter (cold dark matter and baryons) can be neglected and we can describe it as dust 
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with w = 0. Further on, for radiation and neutrinos we have w = 1/3 and finally, for 
the cosmological constant the EOS parameter is = — 1. 

Except for the very early universe, interactions between different species of the 
cosmic fluid do not change their relative number density. Therefore, we can separate 
the density into a sum over independent components, p = J2iPi- Setting the present 
value of the scale factor, oq, to unity (which will be assumed from now on, without 
loss of generality), and evaluating the flrst F-L equation at the present time (t = to); 
we get: 

where Hq (the present-day value of the Hubble parameter) is the Hubble constant. 

We can now deflne dimensionless density parameters of respective components, very 
useful in cosmological analyses: 

n - SttG ^ 
as well as the density parameter of the cosmological constant: 

The quantity pc = (3iJQ)/(87rG) is called the critical density (i.e. the density of a 
flat Universe with matter only, the Einstein- de Sitter model, EiNSTEIN & DE SiTTER 
1932) and the present-day density parameter of non-relativistic matter, Qrnfl = Pm,o/pc, 
will from now on be simply denoted as Q^; similarly for f^A- 

The Universe at low redshifts is dominated by non-relativistic matter: baryons and 
collisionless cold dark matter and, as strongly suggested by observations mentioned at 
the beginning of this Chapter, the cosmological constant (dubbed also dark energy^). 
The relevant equation for the temporal evolution of the scale factor (expansion of the 
Universe) is then 



^In general, the term dark energy has a wider meaning and applies also to generalizations of the 
cosmological constant with an equation of state different than p ~ —p. 
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Figure 1.1: Evolution of the scale factor as a function of the cosmic time, normalized by the 
Hubble constant. Different curves represent different values of density parameters 0,^ and 
as in to (today). [Reproduced from Lahav &: SuTO (2004) by permission of the LRR.] 

This means that the expansion history of the Universe in the Friedman-Lemaitre mo- 
dels can be constrained if we know current values of relevant density parameters^. 
Figure 1.1 shows a(t) for several values of (r2m,f2A). The currently favored ACDM 
cosmological model is the one with fi^ ~ 0.3 and — 0.7 (light-blue dotted line). 



The Friedman-Lemaitre model as outlined in the previous section explicitly assumed 

the homogeneity and isotropy of the Universe. These assumptions certainly hold for 

the early stages of the evolution of the Universe, which is confirmed for instance by 

the high degree of isotropy of the CMB. There are various reasons to believe that this 

■^In the present Universe, both for neutrinos and radiation we have fl^, ~ ~ 0. However, 
Equation (1.10) can be modified to account for their non-zero values. 



1.2 Perturbation theory 
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approximation is valid also currently for the largest scales (e.g. HoGG et al. 2005). 
However, the mere existence of planets, stars or galaxies means that for sufficiently 
small scales, our Universe is strongly inhomogeneous. Cosmology tries to answer the 
question how the structures such as galaxies, their clusters and superclusters, as well 
as voids between them (the 'cosmic web' in general) came to be. The widely accepted 
answer says that they emerged from tiny fluctuations existing at a much earlier epoch 
and grew by the gravitational instability. We will now shortly present this theory, 
especially for small density fluctuations, i.e. in the linear regime. Further details can be 
found in classic textbooks such as Peebles (1980) or reviews, for instance Strauss 
& WiLLICK (1995). 

We will deal with scales much smaller than the so-called Hubble radius, rn = c/Hq, 
and generally with small values of gravitational potential. This means that we can 
focus on the non-relativistic regime, where the Newtonian approximation is valid. For 
a fully relativistic description of the gravitational instability, the interested reader is 
directed to e.g. Mukhanov (2005). 

The basic equations for self-gravitating fluid are given by continuity, Euler and 
Poisson equations: 

^^+v-{pu) = 0, {Tn} 



^ + {u-V)u = --Vp-V$, [132] 
at p 

= AnGp. [lis] 

It is convenient to rewrite the above equations in the so-called comoving coordinates. 
For that purpose we define the position x in the comoving frame, the peculiar velocity 
V, the density contrast 6 and the gravitational potential in the following way: 



X = —— , 1 1.141 

a{t) ^ ' 

V = a{t)x , [lis] 

S{t,x) = ^-1, [He] 

(f){t,x) = ^ + ^aa\x\^, [1.17] 
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where p{t) is the mean density of the background at a given time. In other words, the 
pecuhar velocity i; of a galaxy is the one that adds to the Hubble flow: 

r = dx + ax = H{t)r + v . [l.lS] 

Equations (1.11)-(1.13) in the comoving frame take on the form: 

^+iv.|(i+*H = 0, m 

dv 1 . ^. a 1 „ 1 „ , 

at a a pa a 



72, 



(where the derivatives are now defined in the comoving coordinates). The above equa- 
tions are exact in the Newtonian regime and serve as a starting point for fully-nonlinear 
analyses of gravitational instability. Here we will deal mainly with the linear theory, 
where the density fluctuations are small (5^1). If we now consider only dust and 
neglect pressure, then in the linear regime Eqs. (1.19) and (1.20) can be linearized 
with respect to 5 and v: 

I + = 0, {1^ 

ot a 

ot a a 

Taking now the time derivative of the continuity equation (1.22) and using the Poisson 
(1.21) end linearized Euler (1.23) equations, we get: 

This formula for the linear evolution of the density contrast is a second-order partial 
differential equation in time only, so we can separate temporal and spatial dependencies 
and write the general solution as: 

5 = A{x)Dg{t) + B{x)Dd{t), 

where Dg and Dd are called respectively the growing and decaying mode. The latter 
monotonically decreases as the Universe expands and eventually becomes negligibly 
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small, so we can focus only on the growing one, written from now on simply as D. It 
is given by 

a 7 da' , s 







In the general case of a non- vanishing cosmological constant, this integral is not ana- 
lytic. Closed-form expression for the both modes do exist for A = 0, see e.g. Peebles 
(1980). 

The growing mode can be equivalently expressed in terms of the redshift z: 



z 

In the case of spatially-flat models with the cosmological constant, such as the ACDM, 
it can be written as (e.g. Lahav & SUTO 2004) 



+ 3 / hTf:;3 ) dg, ^ • (128] 





In the regime when the growing mode dominates, we can rewrite Eq. (1.22) as 

V-v = -a(5^ [129) 



and if we define the growth factor f as 



dlnD D 



d m a H D 

we finally obtain the linear-theory relation between the density and peculiar velocity 
divergence fields: 

V -v = -aH fS . (131) 
In the low-redshift universe, 2 -C 1, this relation reads 

V-v = -HofS. [I32] 

The growth factor / is a function of Q^a and A, although the dependence on the 
latter is extremely weak (e.g. Lahav et al. 1991; Martel 1991). A very accurate 
approximation is the one provided by LiNDER (2005): 



/(^]^,A)^fi^(^)^ [133] 
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with 

7 = 0.55 + 0.05(1 +U7|^=i) , [L34l 

where w{z) denotes the equation of state parameter of generahzed dark energy. For the 
cosmological constant, w = —1, so / ~ ^m^^, which reproduces the exact growth factor 
in the ACDM model for any redshift to better than 0.05% over the range G [0.22, 1] 
and it remains accurate to below 1% for > 0.01. A recent paper by BuENO 
Belloso et al. (2011) presents exact solutions in the case of vacuum energy that 
can be parametrized by a constant equation of state parameter w and a very accurate 
approximation for the ansatz w{a) = Wq + Wa{l — a)- 

1.3 Statistical description of cosmological fields 

Galaxies are not distributed in space randomly, but they tend to gather in groups, 
clusters and superclusters, separated by huge voids. This means that the probability 
of finding a galaxy at location x is not independent of whether there is a galaxy in 
the vicinity of x or not: it is more probable to find a galaxy in the neighborhood of 
another one that at an arbitrary location. This is described mathematically in terms 
of the correlation function, ^ (e.g. Peebles 1980). The probability of finding a galaxy 
in a volume element d^ at a location x and at the same time finding another galaxy 
in this volume element at a location y is 

P,, = {ndVf[l+U^,y)] , [L35] 

where h is the average number density of galaxies and C,g{x, y) is their two-point 
correlation function. By analogy, the correlation function for the total matter density 
can be defined as 

{p{x)piy)) = -p'{[l + 5{x)] [1 + 5{y)]) = p ^ [1 + {5ix)5{y))] = p [1 + e(=r, y)] , 

because the mean (expectation) value {6{x)) = for all locations x. In the above and 
hereafter, angular brackets denote averaging over an ensemble of distributions that 
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all have identical properties. Moreover, since the Universe is considered statistically 
homogeneous, ^ can only depend on the difference x — y and not on the individual 
locations. Additionally, owing to the statistical isotropy, E,{x — y) = ^{\x — y\), i.e. it 
depends only on the separation r = \x — y\. Therefore, ^ = ^{r). The galaxy-galaxy 
correlation function and that of matter are in linear theory assumed to be related via 
the linear biasing paradigm: 

Ur) = b''Ur) (137) 

(see e.g. Strauss & Willick 1995), where h is the biasing parameter. We will come 
back to the issue of biasing in Section 1.4. 

An alternative and equivalent description of the statistical properties of matter 
clustering in the Universe is the power spectrum., P{k). It describes the level of structure 
as a function of the length-scale L ~ 27r/k, where is a wavenumber in Fourier space. 
The larger the P{k), the larger the amplitude of the fluctuations on a length-scale 
211 /k. Technically speaking, the density fluctuation fields is decomposed into a sum 
of plane waves of the form 5{x) = J2(^kC0s{x ■ k), where k is the wavevector and 
Ofe is the amplitude. This is a Fourier decomposition of the density field and the 
power spectrum P{k) describes the distribution of amplitudes with equal k = \k\. The 
correlation function and the power spectrum form a Fourier transform pair, that is 

ar) = j^^lp{k)e''^-d'k. (T38] 

In principle, one can also describe the peculiar velocity field in a similar way, 
although this description is more complicated due to the vector nature of velocities. 
This however will not be needed for the purpose of this work. We redirect the interested 
reader to such textbooks as Peebles (1993) or Coles & Lucchin (1995). 

1.4 Cosmic density and velocity fields 

As we have seen in Section 1.2, in the limit of small perturbations, there is a linear 
relation between the velocity divergence and density contrast: 

V•^; = -i/o/(^^m,A)5. [139] 
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We can now invert this Equation via the methods of electrostatics to obtain the integral 
relation, now in proper coordinates: 



An 

(the right-hand side accepts adding an arbitrary divergence-free term, which however 
corresponds to the decaying solution and will be neglected from now on). Writing 
down the peculiar acceleration vector at a position 'r' as 

g{r) = Gp,j 5{r') -^^^ , \1M) 

we immediately obtain the linear relation between (peculiar) velocity and acceleration 
fields: 

= ^^^^ 9W. (mD 

Relations (1.39) and (1.42) can be interpreted in the following way. If we accept 
the gravitational instability as the mechanism of large-scale structure formation, what 
we obtain is that within this framework inhomogeneities in matter distribution in- 
duce gravitational accelerations, which results in galaxies having peculiar velocities 
that add to the Hubble flow. These velocities in turn enhance the growth of the in- 
homogeneities, causing strong coupling between cosmic velocity field and large-scale 
matter distribution. In perturbation theory of Friedman-Lemaitre models, in the lin- 
ear regime, peculiar velocities and accelerations are aligned and proportional to each 
other at every point. More importantly, the proportionality coefficient of this relation 
is a simple function of the cosmological parameter of non-relativistic matter density, 
and practically does not depend on the cosmological constant (nor other forms of 
dark energy). 

However, there are some complications. As what we observe are galaxies, we have 
to assume some relation between their density field and that of matter in general 
(including dark). As we have mentioned in the preceding Section, in linear theory 
this is usually done via the linear biasing paradigm, here expressed in terms of density 
contrasts: 

5^ = h5^. [143] 
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This formula is not strictly valid as for underdensities (voids) with —1 (pm ~ 0) 
and for 6 > 1, we would obtain Sg < —1, which is not possible. The formal definition 
of linear biasing is expressed in terms of galaxy-galaxy and mass-mass correlation 
functions, as was given by Eq. (1.37). The relation (1.43) follows from this, but the 
reverse is not true. This biasing scheme, valid for sufficiently large scales, neglects 
the stochasticity, as well as possible scale- and galaxy-type dependence in the relation 
between the two density fields. For more details on modeling the biasing, see for 
example the review by Lahav & SUTO (2004) and references therein. 

Including the biasing relation into Eq. (1.41) and using the fact that for a spherical 
survey / yZr\:i ~ 0' S^^ following expression for the peculiar acceleration: 

9W = f/ftM^^dV. (144) 

Owing to the nature of astronomical observations, we observe only the distribution 
of luminous galaxies and can measure their peculiar velocities. Using the general 
framework described here, we cannot generally constrain the bias and growth factor 
independently. The two are thus combined into the parameter 



Comparing now Eqs. (1.42) and (1.44), we get the proportionality valid in linear theory: 



v^l3g, (T46] 

or in differential, equivalent, form: 

V -vo^-pS . [147 



These two relations have been widely used to constrain the /3 parameter for more 
than two decades now'^. There are generally two classes of methods for this: the first 



■^The same /3 parameter can be also constrained from so-called redshift space distortions, which 
arise because of the effect of peculiar velocities on the shape of structures seen in redshift surveys. 
A more detailed description is beyond the scope of this work, see e.g. Kaiser (1987) or Hamilton 
(1998). 
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one, velocity-velocity comparisons, is to predict the velocity field from redshift surveys 
using Eq. (1.40) and to compare it with observed peculiar velocities. The second one, 
density- density comparisons, does the opposite: predicts the density field from peculiar 
velocity surveys using Eq. (1.32) and compares it with redshift surveys. For some recent 
results of using the two methods to derive the /3 parameter, see e.g. Erdogdu et al. 
(2006b) or Davis et al. (2011) and references therein. 

Both v-v and d-d comparisons have their advantages and limitations. The latter 
include our poor knowledge of peculiar velocities: they can be measured only if we 
know both the redshift and the actual distance of the galaxy. Apart from the few 
closest galaxies, distances estimated independently of redshifts are constrained with an 
accuracy no better than 10-20%, and only for several thousand galaxies with distances 
smaller than several dozens Mpc. This translates directly to big errors in peculiar 
velocities (bigger at large distances than the estimated velocities themselves). However, 
there is one special peculiar velocity that we know to a very good accuracy of ~5% and 
which can be very well used for cosmological studies. It is the velocity of the system 
that our Galaxy belongs to: the Local Group of galaxies. 

1.5 Local Group as a probe: the clustering dipole 

The Local Group (hereafter 'LG') is a gravitationally bound group of galaxies with 
several dozen members (up to 60, Lee & Lee 2008), including the major players, 
the Milky Way and M31. Apart from internal motions of its galaxies, the whole 
system moves through the Universe with respect to the Hubble flow. Application of 
the Equation (1.41) is simple in this case: if we take r = (the barycenter of the LG) 
and include biasing, we obtain 



This relation could in principle be applied directly to the motion of the LG through the 
Universe. Consequently, comparison of the peculiar velocity and acceleration of the 
LG may serve as a tool to estimate the (3 parameter. Independent knowledge of biasing 
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allows to estimate the cosmological density fim. In reality, however, this procedure is 
not that straightforward. In order to understand why it is so, let us explain how the 
two quantities in this relation, LG velocity and acceleration, are measured. 

The peculiar velocity of the LG is the easier one of the two to constrain. It is known 
from the observed dipole anisotropy of the cosmic microwave background (HiNSHAW 
ET AL., 2009), interpreted as a kinematic effect, and reduced to the barycenter of the 
LG (COURTEAU & VAN DEN Bergh, 1999). It equals to w^mb = 622 ± 35 km/s and 
points in the direction {I, b) = (272° ± 3°, 28° ± 5°) in Galactic coordinates (Hydra 
constellation), where the errors both in amplitude and direction come mostly from 
the uncertainty of the Local Group internal dynamics. The kinematic interpretation 
of the CMB dipole is strongly supported by such properties of the CMB as much 
smaller amplitudes of the quadrupole and higher-order moments of its temperature 
anisotropies, and by the observed alignment of this dipole with the direction of the 
peculiar acceleration of the LG, although the latter is much more difficult to estimate. 
Constraining it requires knowledge of mass distribution in our cosmic neighborhood, 
and its determination had not been possible until deep all-sky galaxy catalogs became 
available. For that reason, the first attempts to measure the acceleration of the LG 
were made not earlier than about three decades ago ago (Yahil et al., 1980; Davis 
& HUCHRA, 1982). 

Using an all-sky catalog, such measurement can be made under the assumption 
that visible (luminous) matter is a good tracer of the underlying density field. The 
general procedure is to estimate the so-called clustering dipole of a galaxy survey and 
infer the acceleration of the LG. However, such inference requires several conditions to 
be met. First, the survey should cover the whole sky; second, the observational proxy 
of the gravitational force (most often the flux of the galaxy in the photometric band of 
the survey) should have controlled properties; and last but not least, the survey should 
be deep enough for the dipole to converge to the final value that we want to find. As 
usually one or more of these assumptions do not hold, the clustering dipole is a biased 
estimator of the acceleration and the estimation of the latter from the former can be 
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done only if the mentioned effects are properly accounted for. 

What is more, in reality we do not observe continuous galaxy density field, but 
instead discrete objects, even if in a very large number. In the following derivation, 
which can be found e.g. in ViLLUMSEN & Strauss (1987), we model galaxies as point 
sources: Pg{r) = X]j S£){r — ri), where 6d is Dirac's delta; Mj and are respectively 
the mass and the position of the i-th galaxy. Putting the coordinate system at r = 0, 
we obtain the acceleration of the LG as a sum of force contributions from all sources 
in the Universe: 

9 = tT.mA, [149] 

(note that as we are interested in the motion of the Local Group as an entire system, 
the galaxies of the LG should not be included in the summation). This Newtonian 
formula^ is still not useful for calculations based on observational data, as masses of 
individual galaxies are known very poorly, or not at all. However, if the i-th galaxy 
has an intrinsic luminosity Lj, we can write 

^~ b b ' ^ 

where Si = Li/Anrf is the flux received from the z-th object. Relation (1.50) means 
that if we know the behavior of the mass-to-light ratio in the band(s) of the survey, we 
can even try to estimate the acceleration of the LG from a two-dimensional catalog, i.e. 
one containing astro- and photometric data only (positions and fluxes). Furthermore, 
if we assume that the mean mass-to-light ratio is a universal constant, T = {M/L), we 
finally get 

i 

In some applications, including the present one, it is more convenient to work in terms 
of matter and luminosity densities. This is especially the case when the luminosity 
density, j, is known for a given band, rather than the mean mass-to-light ratio. We 

^The Newtonian limit can be applied as our whole analysis concerns distances well below the 
Hubble radius rn = c/Hq = 3 Gpc/h. 
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have 



/M\ m^Hl ^— 



which gives 

3 Vtyn Hn 



9 = -^T.s^n. m 

J i 

The luminosity density j for a particular band of the survey can be calculated for 
example from the luminosity function ^{L) of galaxies in this band (e.g. Peebles 
1993) 

oo 

i = j L dL . [1541 



Note that using the Relation (1.53) in Eq. (1.42), we get the linear-theory velocity of 
the LG measured from the flux dipole as 

3 i 

where g represents the acceleration of the LG scaled to units of velocity. The term 
Sifi is the theoretical flux dipole moment of all sources down to the zero flux over the 
whole sky. The universal luminosity density j , measured from a fair sample of galaxies 
in the given band with known apparent luminosities and redshifts, is proportional to 
Ho, which means that the overall result does not depend on the Hubble constant. 

A main complication to the above comes from the fact that realistic galaxy catalogs 
will never reach down to zero flux, irrespective of the used wavelength. Surveys are 
usually flux-limited, and the number of observed sources, A^, is finite. For that reason, 
in the following we will denote the flux dipole of a finite, flux-limited sample as d: 

^^3^Z^ [156) 
2bj Y 

Note that the clustering dipole calculated from the above formula may be a biased 

estimator of the peculiar acceleration of the Local Group, Eq. (1.41). This can be 

overcome by extrapolating the measured dipole to zero flux (Villumsen & Strauss, 

1987). 

The situation becomes more favorable in the case of galaxy redshift surveys. We 
can then use the redshifts as distance estimators and weight galaxies with the inverse 
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of the selection function of the survey (e.g. Yahil et al. 1991), in order to mimic an 
'ideal', volume-limited catalog. However, despite an outstanding advancement in sur- 
veying the cosmos in recent years, the deepest and densest all-sky redshift survey, the 
2MASS Redshift Survey (2MRS, Huchra et al. 2012) contains only ~45,000 galax- 
ies and has a median depth of merely ~ 90 Mpc/h {z-^ ~ 0.03). On the other hand, 
the 'parent' catalog of this survey, namely the Two Micron All Sky Survey (2MASS, 
Skrutskie et al. 2006), reaches 3 times deeper and includes about 20 times more 
galaxies (for details, see Chapter 3). Its redshift coverage, when matched with other 
surveys, such as the Sloan Digital Sky Survey (SDSS, Aihara et al. 2011) or the Six 
Degree Field Galaxy Survey (6dFGS, Jones et al. 2009), is however non-uniform 
both on the sky and in depth (for a recent compilation, dubbed 2M++, see Lavaux 
& Hudson 2011). 

For the purpose of our work, we have thus decided to sacrifice the advantages of 
weighting galaxies, which is possible for redshift surveys, obtaining instead an over- 
whelmingly greater number of sources and unprecedented depth of the survey with 
photometric data only. An additional motivation of using the dipole (1.56), con- 
structed only with the use of fluxes and angular positions of individual galaxies, is 
the fact that it is free of any redshift- space distortions, and in particular of the rocket 
effect (Kaiser, 1987). The latter consists in the fact that the peculiar acceleration 
of the LG calculated using redshifts instead of real distances (which are mostly not 
known) will differ from the actual LG acceleration due to a spurious contribution from 
the galaxies that are in the direction of the LG motion. Here, we do not use distances 
measured in redshift- nor in real space, and the only possible effect of that kind would 
be the anisotropy modulation in the distant galaxy distribution due to the aberration 
effect, which is however completely negligible for our sample (Itoh et al., 2010). 
The only stage at which the Kaiser effect comes into play is in the measurement of the 
luminosity function $(-L) and consequently the luminosity density j. This is addressed 
in the relevant papers where ^{L) is measured, see e.g. Jones et al. (2006). 

The machinery described in this Introduction will be applied in the following Chap- 
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ters of this thesis. We start in Chapter 2 by a case study of a non-hnear effect that 
could possibly be important for our measurements. Subsequent Chapter 3 presents the 
2MASS Extended Source Catalog and how it was prepared for our purposes. Then, in 
Chapters 4 and 5 we apply the data from the 2MASS dataset to calculate the clustering 
dipole of this catalog and show how to use it to measure the (3 parameter. We propose 
two approaches for this: one is to analyze the growth of the clustering dipole with 
increased depth of the sample and compare it with theoretical expectations (Chap- 
ter 4); the other is to use the maximum-likelihood estimation method and optimize 
the window through which the measurement is done (Chapter 5). We conclude and 
present some future prospects in Chapter 6. 
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2 

Non-linear effects: a case study 



^ I "^he perturbation theory described in Chapter 1 apphes to the hnear regime. Math- 
ematically, this means that considered density contrasts must be small, |5| ^ 1 
(as well as should be the velocity divergences). We can thus use the linear theory only 
for sufficiently large scales and for that purpose, non-linear over- and underdensities 
are usually smoothed out in analyses of density and velocity fields. However, struc- 
tures in the Universe nowadays, such as galaxy clusters or big voids, are indeed highly 
non-linear and they exist also in the cosmic neighborhood of the Local Group: the 
Virgo and Coma Clusters, the Shapley Concentration or the Local Void, to name just 
a few. The question then arises in the context of our work: are such objects important 
for the analysis of the clustering dipole as described in Section L5? Can they introduce 
significant systematics? 

There are at least two effects that can be significant here. One is the fact that the 
two vectors v^^^ and g^^^ are not parallel in reality: there is a non-zero misalignment an- 
gle between them. This is both known from observations (e.g. Strauss et AL. 1992; 
SCHMOLDT ET AL. 1999; Maller ET AL. 2003; Erdogdu ET AL. 2006a) and 
expected from simulations (e.g. C1ECIEL4G ET AL. 2001). However, as the observed 
angle is usually small, the deviations from linear-theory relation are usually neglected 
in this respect. There still remains another issue, related to observational constraints 
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that cannot be bypassed: our inability to observe the whole Universe. First of all, 
every galaxy survey is limited (usually by a minimum flux). This causes effects such 
as the shot noise (e.g. Strauss & Willick 1995) related to the fact that only the 
intrinsically brightest galaxies are seen near the edge of the survey. Secondly, irrespec- 
tive of the wavelength in which the observations are performed, our Galaxy obscures 
a significant part of the sky, creating the so-called Zone of Avoidance (ZoA) through 
which we cannot see extragalactic objects. Any 'all-sky' cosmological survey will al- 
ways miss some information in the Galactic plane and its vicinity, due to obscuration 
by dust, gas and stars of the Milky Way. This means that there can always be some 
objects that we are unaware of, but which have an influence on the motion of the LG. 
Should such a hidden object be highly non-linear (i.e. very over- or underdense) and 
close, it could possibly bias the calculated acceleration. 

In order to compute the clustering dipole from the available data, this fact is over- 
come by artificially filling the ZoA, in a more or less sophisticated way (see e.g. Lahav 
1987; Lynden-Bell et al. 1989; Plionis 1989; Maller et al. 2003; Pike & 
Hudson 2005; Erdogdu et al. 2006a; Basilakos & Plionis 2006; Lavaux 
ET AL. 2010). Still, regardless of the method one chooses, such masking of the Galac- 
tic Plane and Bulge will never completely account for any possibly existing, although 
unknown, large-scale structures behind the Milky Way. On the other hand, it is of- 
ten claimed that the direction and amplitude of the calculated clustering dipole do 
not change significantly for different methods of filling the ZoA (e.g. Lahav 1987; 
Erdogdu et al. 2006a). However, this does not necessarily mean that possible 
structures obscured by the Galaxy would have no influence on the calculation of the 
acceleration of the Local Group. 

This issue has been studied by Loeb & Narayan (2008), who analyzed the dy- 
namics of the LG by comparing the dipoles of the CMB and of the 2MRS, focusing 
on the lack of surveyed galaxies behind the ZoA. In order to match the two dipoles, 
they have inferred excess peculiar velocity of the LG towards the ZoA and proposed a 
hidden nearby galaxy or a galaxy cluster as an explanation of the gravitational pull. 
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This suggestion is however in conflict with the current observational knowledge, as 
already in 2004 there was confidence that all significant nearby large-scale structures 
behind the ZoA were known and obscuration of a big galaxy was excluded (Fairall 
& Lahav, 2005). Moreover, LoEB & Narayan (2008) claim that a perfect method 
of filling the ZoA would give 'no discrepancy between the 2MRS and the direction of 
the CMB dipole' (i.e. between the direction of the acceleration and of the velocity of 
the LG). We cannot agree with this statement, since there are other sources of the 
misalignment angle between the two vectors, such as the scatter in the mass-to-light 
ratio (cf. Crook et al. 2007) and stochasticity in the non-linear relation between the 
velocity and acceleration of the LG (Bernardeau et al. 1999; C1ECIEL4G et al. 
2001; Chodorowski et al. 2008). As for the excess pecuhar velocity of the LG, 
TuLLY ET AL. (2008) propose a different explanation (partial at least), basing on 
observational data of distances and velocities of nearby galaxies: motion away from 
the Local Void (LV). 

In this Chapter we will present study of the importance of a non-linear 

and proximate structure, partially hidden behind the ZoA. We choose the LV for this. 
The Chapter is organized as follows. In Section 2.1 we shortly describe the Local 
Void. Section 2.2 presents our calculations of the 'spurious' acceleration induced by 
random filling of the LV behind the ZoA. Next, Section 2.3 covers the issue of the 
directional shift of the clustering dipole due to the discussed effect. Possible corrections 
of the density parameter measured from density-velocity comparisons are addressed 
in Section 2.4. Finally, in Section 2.5 we summarize and conclude. 

2.1 A case study: the Local Void 

We will now analyze the possible influence of the LV on measurements of the clus- 
tering dipole. We would like to emphasize that we do not examine here the importance 
of the Local Void for the very motion of the Local Group. As was shown e.g. by Tully 
ET AL. (2008), the push from the LV is a substantial component of the peculiar ve- 
locity of the LG. In this Chapter we are only interested in the effect of masking the 
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Figure 2.1: Region of the Local Void. The ellipses outline the three apparent sectors of 
the Local Void. The solid dark blue ellipses show two projections of the Inner Local Void, 
bounded on one edge by the Local Sheet. The North and South extensions of the Local 
Void are identified by the light blue short-dashed ellipses and the green long-dashed ellipses, 
respectively. These sectors are separated by bridges of wispy filaments. The red vector 
indicates the direction and amplitude of our motion away from the void. The axes are 
in Supergalactic coordinates in units of velocity. The Local Group is at the origin of the 
coordinate system, at the beginning of the red vector. [Reproduced from Tully et al. 
(2008) by permission of the AAS.] 



intersection of the LV and the ZoA for the purpose of calculation of the clustering 
dipole within linear theory. Our aim is to investigate possible systematics related to 
the fact that a part of the LV is hidden behind the ZoA. This will be useful in assessing 
the total errors of the measurements presented in Chapters 4 and 5. 

The Local Void is a structure first identified in the Nearby Galaxies Atlas (Tully 
& Fisher, 1987). Although our close cosmic neighbor, it is poorly defined because 
much of it lies behind the plane of the Milky Way: its center is located at Galactic 
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coordinates / ~ 30°, 6^0° (Tully 2007, private communication; but see Nasonova & 
Karachentsev 2011 and references therein for other determinations). It is neverthe- 
less confirmed (by surveys in different wavelengths) that there is an under abundance, 
though not a total lack, of galaxies in a very large part of the sky at low redshifts 
(e.g. Kraan-Korteweg et al. 2008). This empty region begins at the edge of the 
Local Group, with the so-called Local Sheet as a bounding surface (Tully 2008a, B; 
Tully et al. 2008). Both the shape of the LV and its dimensions are not exactly 
known; however, it seems to have two components: a smaller void with a long dimen- 
sion of ~ 35 Mpc enclosed within a larger void with a long dimension of the order 
60 ~ 70 Mpc (Tully 2007; Tully et al. 2008). Figure 2.1 (reproduced from 
Tully et al. 2008) presents the observational data in the region. 

In the framework of cosmological gravitational instability, voids are expected to 
form out of initial under densities, i.e. regions less dense than average (e.g. Hoffman 
& Shaham 1982). Their expansion is faster than the Hubble flow, which results in 
voids 'swelling'. Tully et al. (2008) calculated an 'effective Hubble rate' of the 
Local Void, under the simplest assumption of it being empty and spherical, and used 
this calculation, together with the peculiar velocity of the LG away from the LV, to 
estimate the radius of the latter as at least 16 Mpc. However, these dynamical esti- 
mates cannot exclude an effective diameter of the LV greater even than 45 Mpc (Tully 
2008a,b; Tully et al. 2008). 

Such a prominent structure, partly hidden behind the Zone of Avoidance, should 
have an influence on the calculation of the integral in Eq. (1.48) and consequently, on 
determinations of /3 and from the - S'lg relation. In the following Sections, 
using a simple model, we will analyze the amount of the systematic error it may cause 
for the estimation of the acceleration of the Local Group. 
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2.2 Acceleration due to (the lack of) the Local Void 

Let us define the scaled acceleration vector g^^ in convenient units as 

(we omit at the moment the normalization, not important in the qualitative analysis 
that follows). We can split the integral into two parts: one over the Local Void and 
the other covering 'the rest of the Universe': 




LV R3_Lv 



For simplicity, we write only 'LV in this and the following formulae, although we mean 
in fact 'LVnZoA', as will be explained later in detail. Now, if we make a simplifying 
and maximizing assumption that the Local Void is completely empty (plv = 0), we 
have (5lv = —1- On the other hand, if we filled the part of the Local Void hidden 
behind the ZoA by randomly chosen 'galaxies', i.e. if we assumed average background 
density in this part (plv = Ph), we would get 5lv = 0. We would thus measure some 
acceleration that we will call spurious, which is the difference between the calculation 
made with the ZoA filled randomly and the true peculiar acceleration of the Local 
Group: 

^spur = ^cal — St = (2.3] 
LV R3-LV LV ' ' R3-LV LV ' ' 

Let us now calculate explicitly the above integral, returning to physical units. 

2.2.1 Spherical void 

For that purpose, in addition to the assumption of the Local Void being completely 
empty, we start by modeling it as a sphere with some radius R, and assume that the 
Local Group lies exactly at the edge of the LV (see Figure 2.2)^. We place the center of 
^We neglect the size of the Local Group, estimated usually to about 1.5 Mpc (van den Bergh, 

2007; VAN DER MAREL & GUHATHAKURTA, 2008). 
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Figure 2.2: A schematic view of the Local Void: (left) along the z axis, with the Galactic 
plane face-on; (right) along the y axis, with the Galactic plane edge-on. On both panels, the 
relative sizes of the Local Group and the Local Void are not to scale. 



the Local Void in the plane of the Galactic equator, in accordance with observations, 
which additionally maximizes the studied effect. The x and y axes of the coordinate 
system lie in the Galactic plane, with the origin in the center of the LV (Figure 2.2, left 
panel). As the latter is situated at /lv = 30°, 6lv = 0°, this means that our coordinate 
system is shifted with respect to the Galactic one and rotated by 210° in the Milky 
Way plane (our z axis is parallel to ^cai)- The part missed in all-sky surveys due 
to the Zone of Avoidance is a spherical section bounded by two planes (Figure 2.2, 
right panel); the angle between them, a, depends on the survey, but for our purposes 
a ~ 20° (we consider near- infrared wavelengths, as is explained hereafter). Owing to 
the smallness of this angle [a ~ 0.35 rad), we treat the masked region as a thin wedge. 

We additionally include a shell of compensation at the edge of the LV, of width 
AR. According to the standard picture of void formation, where such structures are 
grown from initial underdensities and gradually expand faster than the cosmological 
background, the matter expelled from inside of the void creates a layer at the edge, 
of density higher than the average background one (e.g. Maeda et al. 2011 and 
references therein). This is supported by both simulations (see e.g. van de Weygaert 
& ScHAAP 2009) and observations (cf. Fairall 1998 and references therein): voids 
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are usually bounded by 'filaments' and 'walls' of high density contrast. In our cosmic 
neighborhood, a part of such a shell is probably the Local Sheet and the Local Group 
is located at its inner edge (Tully 2010, private communication). 
A mass element dM contributing to the spurious acceleration is 

dM = p^dV = p^hdS = pi^ardS , (Z4j 

where ~ ar is the height of the wedge at a distance r from the LG and dS is a 
surface element. The differential acceleration 'measured' due to random filling of the 
ZoA will be 

dM{r)r p^^ardSr qdqd(p{R—q) 
dg = — G r = — G = —a G pb 



J.3 j,3 |i? — 

with dS = qdqd(p, q = {q cos (p,q sirup), R = (i?, 0) and < (p < 2n is the azimuthal 
angle. It is easy to check that from symmetry Qy = 0. The x-component of the spurious 
acceleration is given by 

g, = -aGpkRI, {tI} 

where 

J J t + 1 — 2t cos 09 

Thus the value of the peculiar velocity caused by the spurious acceleration is (in linear 
theory)^ 



47rGpb ^ 4' 

It is interesting to compare this value with the one induced by a sphere of radius R and 
density pb- The peculiar acceleration at the surface of the sphere is then g, = ^iiGphR, 
which gives the linear peculiar velocity of 

v. = ^Hof{Qr.)R. 



We thus have 



^spur 3 



V, 4 



- a 



^Herc and later in this Chapter, we neglect the biasing (cf. Eq. 1.43), or in other words, we set it 
equal to unity. Its exact value does not affect the conclusions of this Section. 
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which shows that the effect is expected to be small. 

The above calculation applied to an isolated void without any shell of compensation. 
In order to verify the effect of such a shell, let us first assume that the average matter 
density inside the layer is constant, psh > Pb, and its thickness is AR. The Local 
Group, of negligible size, is still located at the edge of the Local Void, as in Figure 2.2. 
Putting now the mass element inside the intersection of the shell with the ZoA, we can 
see that the differential spurious acceleration will read 

dc/sh = aG (p,h - Pb) \R-q\'i ' i 2-llJ 

with q and R = {R, 0) defined as before, but now R < q < R + AR. The ^/-component 
of the acceleration vanishes from symmetry, and along the x direction we have 



9x 



(sh) 



a G {psh - ph) R^sh , [2.12] 



with Xsh different from X only by the limits of integration in C, and equal to 

^ I AR 

J J t'^ + \ — It cos 09 

1 7- 

This result means also that the spurious acceleration of the LG due to the part of the 
compensating shell hidden behind the ZoA will equally vanish if the density distribution 
Psh is not constant, but depends on the distance from the center of the LV (as we can 
divide the shell into infinitesimally thin layers of constant density each). The bottom- 
line is that a compensating shell with radial density distribution will not affect the 
spurious acceleration of the Local Group provided that the latter is located at the 
inner edge of the shell, as observational constraints suggest.'^ On the other hand, if 
the LG was placed in the interior of the compensating layer, or at its outer edge, the 
discussed spurious acceleration from the LV would be diminished., by up to 2/3 in the 
limit of an infinitesimally thin shell. This means that the calculations presented so far 
concern a maximizing case. 



■^In fact, the spurious acceleration from the hidden part of the shell would vanish even if the LG 
was inside the LV, as long as the center of the latter was coplanar with the Galactic equator. 
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Having calculated the amplitude of the spurious acceleration in the model, we 
can apply our results to observational data. As an example, we take the analysis of 
Maller et al. (2003) of the clustering dipole of the Two Micron All Sky Survey^, 
since these are the same data that will be used in our analysis of Chapters 3-5. Maller 
ET AL. (2003) used the measurements in the near-infrared band and defined the 
Zone of Avoidance as the region with \h\ < 7° for / > 230° or / < 130° and \h\ < 12° 
for / > 330° or / < 30°. Let us take a 'mean value': 

a = 2A6 = ^(2x7° + 2x 12°) = 19° ~ 0.332 rad. [214] 

We now use Eq. (2.8) and for consistency with Tully et al. (2008), we apply 
= 74 km/s/Mpc, = 0.24 and R = 16Mpc, to obtain an approximate value of 
the additional, spurious velocity of the Local Group, measured if the Local Void is not 
properly accounted for: 

^^spur - 45 km/s . [2.15] 

When compared to the velocity of the LG relative to the CMB reference frame, Ucmb ~ 
622 ± 35 km/s, one can see that this effect is of the same order as the error in the 
measurement of v^^^ . Note also that due to the proportionality of Vspnr to the radius 
of the LV (Eq. 2.8), increasing its diameter to 45Mpc, while preserving sphericity, 
would cause the spurious velocity to raise significantly to ~60 km/s. However, owing 
to the geometry of the problem, this would not largely affect general conclusions of 
our analysis (see Sees. 2.3 and 2.4). Moreover, observational constraints point rather 
to some degree of elongation of the LV than to a larger size of the whole structure. We 
will now address this issue. 

2.2.2 Elongated void 

The calculations so far assumed a simplistic model of a spherical void. However, we 
can see for example in Figure 2.1 that the Local Void should be possibly modeled by a 
"^We will come back to the issue of the 2MASS clustering dipole when presenting our own mea- 
surements (Chapters 4 and 5). Here however we want to obtain an estimate of systematics that can 
be used for our calculations, hence we use externally provided earlier results. 
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more sophisticated structure, such as an elhpsoid. Current observational data suggest 
that the LV is elongated in the Supergalactic SGY = plane, roughly coincident with 
the Galactic plane and the x axis of our coordinate system. One should bear in mind 
however that this effect could be a manifestation of the existence of the ZoA itself: 
lack of observed galaxies in this region of the sky may be simply due to obscuration. 
Nevertheless, in case the elongation is real, for completeness of our analysis let us 
examine this possibility. 

For that purpose we assume that the section of the LV in our xy plane is an ellipse 
with semiaxes B = R and A = kR, where elongation k > 1. The major axis of the 
ellipse is placed along the x axis of the coordinate system. As in the spherical case, 
owing to symmetries of the problem, the only non- vanishing component of the spurious 
acceleration 'acting' on the LG is the x one. It is easy to check that the relevant formula 
for Qx is now 

g^c = -a G py, K RIe , [2.16] 

with the integral Xg given by 

J J f + 1 — 2f cos 09 

7- 

Here, = 1/ (^k ^/l — cos^ ip^ is an appropriately normalized equation of the 

ellipse in polar coordinates with the eccentricity e = vT^^/t^. 

The spurious velocity of the LG due to such elongation of the LV will be larger in 
comparison to the spherical case by a factor 

spur -'-e -'-e 
T = , , , = K = K . 



VspUT 



(sph) -T „ 



Obviously, for k = 1, we have t = 1. What is important here is that a linear increase 
in K results only in a slower-than-linear raise of the r factor: for example k, = 2 gives 
r = 4/3 and if k = 3, then r = 3/2. This means that for the observationally allowed 
elongation of k ^ 2 and minor semiaxis i? ~ 15 Mpc (TULLY ET AL., 2008), the 
spurious velocity of the LV would raise by 1/3 to ~60 km/s, which is the same as in 
the previously discussed case of the enlargement of a spherical and empty void. Note 
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however that observations clearly show that the Local Void is not completely empty 
(cf. Figure 2.1) and this high value of the spurious velocity will be an upper limit for 
our considerations. 

2.3 Shift of the clustering dipole 

Knowing the amplitude of the spurious velocity induced by randomly-filled inter- 
section of the Local Void and the Zone of Avoidance, we would like to check the shift of 
the direction of the measured clustering dipole when the effect of the LV is accounted 
for. From Eq. (2.3), the true velocity of the LG (proportional to its acceleration in 
linear theory) is related to the calculated ('measured') one and the spurious component 
via 

■"T = -ycal - '"spur • [2.19] 

The vector i^spur is directed to the center of the Local Void;^ as fx we adopt v^y^^ = 
622 km/s. For the measured clustering dipole we choose the dipole of the 2MASS 
survey, as given in Maller et al. (2003) for random filling of the ZoA: /cai = 266°, 
&cai = 47°. Using these values altogether, after some calculations we find that the 
direction of the dipole is shifted down by 5° in I and 2° in b. Thus, the 'true' direction 
of the 2MASS dipole would be 

It ^ 261° , 6t ^ 45° . [Z20) 

The shift is small, but as a systematic effect, it should be in principle accounted for 

in the measurement of the clustering dipole. However, random filling is not the only 

possible, nor the most optimal, way to deal with the ZoA. A better method, preferred 

by us (see Chapter 3) is to clone the sky below and above the ZoA, which has the 

advantage of approximately tracing structures through the ZoA. For the latter method, 

Maller et al. (2003) obtained /cai = 263°, 6cai = 40°, so differences between the two 

^Note that although in our coordinate system the only non- vanishing component of the spurious 
velocity is the x one, when projected on Galactic coordinates it has xcai and ycai components of 
comparable values, equal respectively to 0.87 Wspm- and O.Swgpur- 
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methods give A/^ai = 3°, Afe^ai = 7°. Therefore, the differences in the direction of the 
2MASS dipole resulting from distinct methods of treating the ZoA are comparable to, 
or even greater than, the effect of the LV. These conclusions do not change significantly 
even if we include the elongation of the LV: for t^spur = 60 km/s, we have a shift by 
Al = 6° and Ab = 2° towards h = 260°, 6t = 45°. 

The shift of the direction of the 2MASS clustering dipole changes the misalignment 
angle with respect to the CMB dipole. However, the amplitude of the shift is compa- 
rable to the uncertainty of the CMB dipole direction (3° and 5° respectively for / and 
b). Moreover, owing to the specific geometry of the problem, the calculated change of 
the misalignment angle turns out to be smaller than 1° even for high (but reasonable) 
values of fspur- We thus conclude that masking the LV has negligible impact on the 
misalignment angle between the 2MASS and CMB dipoles. 



2.4 Correcting the density parameter measurement 

Application of Equation (1.48) serves as a method to measure the cosmological 
parameter fl^, in principle by comparing the velocity of the LG (equal to '^cmb) 
gravitational acceleration inferred from an all-sky galaxy survey (but see the discussion 
in Section 1.5). From Eq. (1.48) it follows that v^^^ = Q^^^ g^^ (with b = 1, cf. 
footnote 2), where g^^ is the gravitational acceleration of the LG expressed in units 
of velocity. The acceleration measured without the LV accounted for results in the 
'calculated' value of Qrn, such that ^'^fi^ = fcMs/^cai- If the LV is taken into account, 
we find the 'true' value of Q^n, i-e. ^t^^ = fcMs/^T- When we divide fi^^f^ by Q^^^, all 
the scaling factors relating the velocity to acceleration in linear theory (including the 
biasing parameter) cancel out. Therefore, 

00.55 

"cal 



1^0.55 



{221} 



l"f 'Veal 

The velocity vt is related to Wcai by Equation (2.19), where Vs^nr is a small correction. 
Therefore, we can expect that the relative change in the value of fim will be small. 
We thus write f2cai = + Afi and expand the expression VP^-^^ l^i:^^ ^o first order in 
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AQ/Q. The result is 



1 + 0.55—. [^22] 



The right-hand-side of Eq. (2.21), calculated using Formula (2.19) with t^spur = 45 km/s, 
is 1.029. Hence finally 

^ = 0.053. [Z23] 

In other words, not accounting for the existence of the LV in measurements of the 
clustering dipole biases the estimated value of by about 5% for the radius of the 
spherical LV equal to 16 Mpc. If we allow for non-sphericity, this bias rises to some 
7%. This means that the infiuence of the LV (as discussed here) on the determination 
of the cosmic density parameter from the comparison of the velocity and acceleration 
is small. Indeed, typically the uncertainty of the degenerate combination /3 = /(f2m)/6 
(with h not necessarily unity), amounts to at least 10~20% (for determinations using 
2MASS data see e.g. Pike & Hudson 2005; Erdogdu et al. 2006a or Davis 
ET AL. 2011). Owing to our ignorance of the exact value of 6, which typically has 
errors as big as 20% (Maller et al., 2005), we can conclude that the total relative 
error in fim is likely to be higher than in /3 and the inclusion of the effect of the Local 
Void would not contribute largely to the final error budget, although it possibly should 
be included systematic effect. 

The above analysis, although concerning a non-linear effect, was performed within 
the linear theory, in which the estimated peculiar velocity, compared to the observed 
one, is inferred directly from the scaled peculiar acceleration (the clustering dipole). 
However, a completely empty void is a non-linear structure (e.g. BiLlCKl & Chodor- 
OWSKI 2008) and the actual spurious velocity of the LG, generated by the LV, will 
be greater than the corresponding scaled spurious acceleration. It is known that for a 
void with (5 = — 1 (i.e. p = 0), the relation is fgpur — l-5(7spur (Bernardeau et al. 
1999; BiLlCKl & Chodorowski 2008). Application of this result would increase the 
systematic effects considered in this Chapter by roughly 50% and enhance their signif- 
icance. Nevertheless, in order to make this approach self-consistent, one would have 
to take into account non-linear effects from all other sources, especially those nearby. 
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both over- and underdense. This would be a very difficult task to perform analytically, 
if not impossible, and is beyond the scope of this analysis, which deals with a simple 
model of the Local Void surrounded only by a shell of compensation. 

2.5 Summary and conclusions 

A serious problem plaguing determinations of the clustering dipole is that every 
survey called 'all-sky' misses a significant amount of galaxies due to obscuration by 
dust, gas and stars in the disk and bulge of the Milky Way (the Zone of Avoidance). 
To overcome this problem, in order to calculate the clustering dipole of the given 
survey, the ZoA is filled with mock galaxies. Their properties are chosen in a way to 
mimic the true, although unknown, galaxy distribution in the obscured part of the sky, 
extrapolating the one known from the rest of the celestial sphere. 

A part of the ZoA intersects with a nearby void region, the Local Void. When the 
existence of such a structure is not accounted for in the calculation of the acceleration 
of the LG, a spurious term is generated. In this Chapter we have calculated both the 
amplitude and the direction of this spurious acceleration. For simplicity, we have ffist 
assumed that the LV is spherical and for its size we have adopted the value estimated by 
TuLLY ET AL. (2008). We have also made the assumption that the LV is completely 
empty. Even then the amplitude of the spurious component amounts only to 45 km/s 
in units of velocity. Including the observed elongation of the LV enhances this value by 
1/3. On the other hand, possible presence of massive structures inside the LV, hidden 
behind the ZoA, could only lower this value. 

This artificial acceleration changes also the direction of the calculated clustering 
dipole. We have shown that this change is comparable to the uncertainty in the 
direction of the peculiar velocity of the LG, determined from the dipole component of 
the CMB temperature distribution, reduced to the barycenter of the LG. Moreover, by 
chance it points almost perpendicularly to the misalignment vector (i.e. the difference 
between the vectors of the velocity and acceleration of the LG). This results in a 
negligible shift of the misalignment angle, by less than 1°. 
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The final effect that we considered was the error in the inferred value of the non- 
relativistic matter density resulting from the negligence of the LV. We have esti- 
mated the relative error in this parameter due to this effect as approximately 5%~7%. 
Therefore, up to this accuracy the influence of the Local Void on the determination 
of from velocity-density comparisons can be neglected. On the other hand, this 
additional systematic should be taken into account in the total error budget of the 
density parameter determined by such a method. 

We would like to reiterate that our results do not negate the dynamical influence 
of the Local Void on the Local Group; on the contrary, Tully et al. (2008) have 
shown that this influence is significant. It is only the effect of masking the intersec- 
tion of the LV and the ZoA that seems to be of little importance for the purpose of 
calculation of the clustering dipole within the linear theory. This partially supports 
the claims that the Zone of Avoidance is not a crucial issue in determinations of the 
peculiar acceleration of the LG from all-sky surveys, especially such as 2MASS, where 
Galactic extinction is much weaker than in optical wavelengths. We will include the 
effect discussed here as a systematic in the error budget of the 2MASS clustering dipole 
analysis, presented in subsequent Chapters. This analysis will start with a presentation 
of the data that we used. 
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The data: Two Micron All Sky Survey 



he analysis of the clustering dipole requires a catalog of galaxies that covers 



the whole sky, and desirably has uniform properties, such as photometry and 
astrometry. These are generally demanding requirements and only several catalogs 
of that type have been obtained so far. Apart from the attempts to gather all-sky 
galaxy data by compiling several sources (e.g. Lahav 1987 or HUDSON 1993), two 
approaches are possible. The first one is to use a single instrument in space; the second 
is to have two identical telescopes on the two Earth hemispheres. The former approach 
was taken by the Infrared Astronomical Satellite (IRAS) team, which resulted in the 
first really uniform all-sky galaxy catalogs (e.g. FiSHER ET AL. 1995), widely used 
afterwards for cosmological measurements, including those of the clustering dipole and 
LG acceleration (see the list of references in the beginning of Chapter 4). The second 
method to obtain extragalactic data for the whole sky, which consists in using two 
ground-based instruments on both hemispheres, was the key idea behind the Two 
Micron All Sky Survey. 



(2MASS) 
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3.1 2MASS Extended Source Catalog 

The Two Micron All Sky Survey (2MASS, Skrutskie ET AL. 2006) is the first 
near-infrared survey of the whole sky (covering 99.998% of the celestial sphere), and was 
performed in the period 1997-2001 in the J(1.25/im), H {1.65fim) and Kg (2.16/im) 
bands, with the use of twin 1.3-m ground-based telescopes, one at Mount Hopkins 
(Arizona, USA) and the other at Cerro Tololo Inter- American Observatory (Chile). 
All the data from the survey are available through the NASA/IPAC Infrared Science 
Archive.^ The main outcome of this project are two photometric catalogs: of point 
sources (PSC), containing about 471 million objects (mainly stars and some quasars), 
and of extended ones, with more than 1.6 million objects, majority of which are galaxies 
(> 98%) with additionally some diffuse Galactic sources (Jarrett, 2004). The 
Extended Source Catalog (XSC), which was used for the purpose of our analysis, is 
complete for sources brighter than Ks — 13.6 mag (~2.7mJy) and resolved diameters 
larger than ~ 10 - 15". The near-infrared flux is particularly useful for the purpose 
of large-scale structure studies as it samples the old stellar population, and hence 
the bulk of stellar mass, and it is minimally affected by dust in the Galactic plane 
(Jarrett, 2004). An additional advantage of using 2MASS data, especially in the 
context of calculating the flux dipole (Eq. 1.56), for which apparent magnitudes are 
used, is the global photometric uniformity of the catalog, which was enforced by nightly 
photometric calibration to an extensive set of standard star fields. This procedure 
assured equally that there is no bias or offset between the photometry or astrometry 
obtained with the two telescopes used for the survey (Skrutskie et al., 2006). 
On the other hand, as any survey, 2MASS is not perfect. It is biased against optically 
blue and low surface brightness galaxies, such as dwarfs, but sensitive to the early type, 
bulge-dominated ones. However, as the former have very small luminosities and masses, 
their possible underrepresentation in the catalog should not influence significantly our 
results. What is more, these biases are much less severe than it was in the case of IRAS 

^http : //irsa . ipac . caltech . edu/Missions/2mass . html 
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which, due to longer wavelengths (middle and far infrared), was sensitive mainly to 
starburst galaxies. 

Figure 3.1, courtesy of Thomas Jarrett^, presents an Aitoff (equal-area) projection 
of the 2MASS XSC galaxy distribution in Galactic coordinates, color-coded by redshifts 
or their estimates from Jarrett (2004). Several local large-scale structures are 
clearly visible, including the Supergalactic plane (right of the 0-th meridian). The 
Milky Way galaxy, creating the Zone of Avoidance superimposed on the 'cosmic web', 
covers only a small fraction of the sky when compared to surveys in other bands. 
Note that the sample gathered by Jarrett (2004), which is a promising first step 
towards deriving photometric redshifts for all those galaxies in the 2MASS XSC that 
do not have spectroscopic ones measured, is not useful for the purposes of our work. 
Those redshift estimates are only preliminary, do not have controlled errors and their 
accuracy depends strongly on the position of the galaxy on the sky. 

3.2 Data preparation 

The 2MASS photometry offers several types of 'magnitudes' for extended objects, 
depending for instance on the type of aperture used. Throughout the whole analysis we 
use the 20mag/sq." isophotal fiducial elliptical aperture magnitudes, which are defined 
as magnitudes inside the elliptical isophote corresponding to a surface brightness of 
yUband = 20 mag/ sq.". We prefer those to the Kron ones as the latter use large and noisy 
apertures, prone to contamination, resulting in systematic overestimation. Our choice 
is additionally supported by the considerations in the appendix of KOCHANEK ET AL. 
(2001). However, we must remember to correct the values used by adding an offset of 
A = —0.2 mag when converting to flux, in order to compensate for the flux lost outside 
the aperture (typically ~ 10% - 20%, Jarrett et al. 2003). We have checked that 
this offset is roughly equal to the one between isophotal fiducial elliptical aperture 
magnitudes and the 'total' ones, obtained from fit extrapolation (the latter are also 

^http : //web . ipac . caltech. edu/staf f /jarrett/lss/ index.html 
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Figure 3.1: Distribution of galaxies from the 2MASS Extended Source Catalog. 
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noisy and not recommended for such analyses as ours). The magnitude correction by 
a constant factor certainly introduces some scatter in total flux estimates, as it may 
depend on galaxy morphology. The latter is very hard to constrain from 2MASS data, 
we will thus treat this scatter as a systematic effect that needs to be included in the 
error budget. Our tests have shown that this error is of the order of a couple percent. 

In order to prepare the data for our purposes, we have proceeded as follows. First of 
all, we applied the extinction correction from Schlegel et al. (1998), by calling the 
procedure dust_getval'^ for Galactic coordinates of each of the objects. The procedure 
yielded values of E{B — V), which were subtracted from the original magnitudes with 
appropriate multiplicative factors Ry taken from Cardelli et al. (1989): 0.902 for 
J, 0.576 for H and 0.367 for K. We performed the subtraction for objects with \b\ > 5° 
due to the statement of ScHLEGEL et al. (1998) that for \b\ < 5° the predicted red- 
denings should not be trusted (the sources in the ZoA were not included in our catalog 
apart from several brightest; see below). Moreover, for some minor parts of the sky, the 
extinction correction gave unreliably high reddenings, which resulted in some objects 
becoming unrealistically bright (with negative magnitudes) and eventually deleted. At 
this stage, we have also identified and removed the following sources (with some found 
in more than one category): 

• artifacts: flag cc_f lg=a in the 2MASS XSC (122 objects); 

• sources with NULL or unreliable K magnitudes (as described above) (718 objects); 

• non-extended sources: flag vc=2 in the 2MASS XSC (7383 objects); 

• Local Group galaxies, taken from the list of Lee & Lee (2008) (31 objects);^ 

• some of the Milky Way sources, taken from a list of 4454 such objects, separately 
identified earlier in the 2MASS XSC by Thomas Jarrett (private communication). 

^http : //www . astro .princeton . edu/~schlegel/dust/dustpub/CodeC/README. C 

''Not all the objects from the Local Group were found in the database. These were some dwarf 
galaxies of low mass and ncar-IR luminosity, hidden behind the Galaxy or with surface brightness 
below the threshold of 2MASS. 



39 



CHAPTER 3. THE DATA: TWO MICRON ALL SKY SURVEY (2MASS) 




Figure 3.2: Histogram showing the J — K color distribution of the 2MRS galaxies (red filled 
bars) and of the 2MASS extended objects with \b\ > 15° (blue open bars). The counts are 
normalized to a common value for clarity. 



3.3 Removal of Milky Way sources 

The 2MASS Extended Source Catalog contains mainly galaxies; however, it is also 
comprised of Milky Way entities, such as stellar clusters, planetary nebulae, HII re- 
gions, young stellar objects and so on. In order to keep our analysis reliable, these 
objects had to be removed from the catalog. This was partially done for the 4454 
sources mentioned above. However, owing to the size of the catalog, any further 'man- 
ual' procedure of Galactic object removal was impossible and only a method based 
on some general properties could be applied. A useful one in this regard is the color, 
i.e. difference of magnitudes in two bands. In their analysis, Maller et al. (2003) 
made a cross-correlation with galaxies spectroscopically confirmed from the Sloan Dig- 
ital Sky Survey (SDSS) and excluded extended sources brighter than K = 12 mag 
with colors J — K < 0.75 or J — K > 1.4; at fainter magnitudes only those objects 
with J — K < 0.5 were removed. We have checked these conditions by examining the 
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J — K distribution of galaxies in the 2MASS Redshift Survey (which by construction 
contains only extragalactic sources) . We have found that indeed galaxies are clustered 
around J — K = 1 (cf. Fig. 3.2); however, we have decided to alter the limits given by 
Maller et al. (2003). Analyzing additionally the distribution of XSC objects with 
K < 13.5 mag and \h\ > 15°, among which there are mainly non-Milky Way sources, 
apart from some molecular clouds (private communication of Janusz Kaluzny), we 
have decided to keep in our catalog those objects that have 0.6 < J — K < 2.0 (see 
also Jarrett 2000). Figure 3.2, showing the histogram of the J — K color for the 
2MRS sample and for the above mentioned 2MASS objects off the Galactic plane, 
supports our choice. We use our criterion for all sources, independently of magnitude, 
as we think that a differentiation as in Maller et al. (2003) could lead to a bias in 
the sample. An additional visual verification of the 100 brightest objects which pass 
this filter off the Galactic plane confirms that indeed all of them are galaxies and that 
only one galaxy with an extreme value of J — is removed by this procedure up to 
Ks ~7.5 mag. 



3.4 Zone of Avoidance 

As was already discussed in Chapter 2, an important issue in the calculation of the 
clustering dipole is the Zone of Avoidance (ZoA), i.e. the region of the sky with small 
Galactic latitudes b, which obscures galaxies behind the Galactic plane and bulge. 
Although the Galactic extinction is much lower in the near infrared than in visible 
bands (Cardelli et al., 1989) and this applies equally to the ZoA (Jarrett 
ET AL., 2000), the 2MASS XSC is still incomplete near the Galactic equator, mainly 
due to high stellar density in this region of the sky (Kraan-Korteweg & Jar- 
rett, 2005). For that reason, and owing to inapplicability of the extinction maps of 
SCHLEGEL ET AL. (1998) for |6| < 5°, we have masked out the Galactic plane and 
bulge in the following way. For the shape of the mask we have chosen the one proposed 
in Erdogdu et al. (2006a), i.e. we have skipped all the objects with |6| < 5° (plane) 
and \b\ < 10° for / < 30° or / > 330° (bulge). Then we have filled the resultant gap 
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Figure 3.3: Schematic illustration of masking the Zone of Avoidance and cloning the adjacent 
strips. Galaxies from the navy areas were removed and replaced by those from above and 
below, with their latitudes mirrored using the edges of the ZoA as the axes of symmetry. 
The observer is located at the center of the circle. 

by cloning the adjacent strips, with mirror-like reflections: for instance, objects with 
10° < b < 20° were copied to the bulge by assigning 6ncw := 20° — 6oid and keeping 
other parameters unchanged (such as the longitude / and magnitudes). An analogous 
procedure was used for the negative latitudes and for the Galactic plane. Figure 3.3 
illustrates schematically this procedure. Such cloning has the advantage over random 
fllling (considered both in Maller et al. 2003 and Erdogdu et al. 2006a and 
discussed in Chapter 2) that it extends the structures from above and below the ZoA; 
moreover, the only artiflcial discontinuity of the galaxy distribution created in this 
procedure is at the Galactic equator and at the edges of the box masking the bulge. 
We have also tried other masks and methods of fllling the ZoA, and found no special 
importance for the results of the analysis presented here. This will be shortly addressed 
in Section 4.1. 
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3.5 Large Galaxy Atlas 

Once the ZoA has been masked and filled, we have added to our catalog several 
galaxies that were not present in the 2MASS XSC but could be found in the 2MASS 
Large Galaxy Atlas (LGA, Jarrett et al. 2003). This atlas^ contains the ~ 600 
largest galaxies as seen in the near-infrared, of which around 50 are not present in the 
XSC or are located in the ZoA, |6| < 5° (17 sources). Among the latter, three are 
of particular importance for the Local Group motion, namely Maffei 1, Maffei 2 and 
Circinus. We will discuss their infiuence on our results later in the text. Note that this 
addition of LGA galaxies does not spoil the photometric uniformity of the resulting 
sample because all the galaxies from the LGA present also in the XSC were assigned 
the magnitudes from the former catalog when the final version of the latter one was 
constructed. 

In case of those LGA galaxies that were present in the ZoA, the Schlegel et al. 
(1998) maps are known to overestimate the Galactic extinction by roughly 15% (e.g. 
Schroder et al. 2007); we have thus decreased the E{B — V) by that amount there. 
The exceptions are Maffei 1 and Maffei 2, for which we used exact values of extinction, 
given in FiNGERHUT et al. (2007), as well as Circinus with E{B - V) = 0.677 
(For et al., in preparation). Note however that apart from those three galaxies, which 
have extinction-corrected magnitudes below K = 5 mag, all the remaining ones added 
from LGA are much fainter, by 2 mag or more, and possible misestimation of their 
extinction does not largely influence our analysis. 

3.6 The final catalog 

The flnal catalog contained 1,464,028 galaxies up to the 14th magnitude in the 
Ks band. Among these, 108,023 were clones copied to the ZoA. Appendix A provides 
information about the 50 brightest galaxies of the sample (excluding the Local Group). 

^Accessible through IPAC Infrared Science Archive (IRSA) at 
http : //irsa. ipac . caltech. edu/applicatioiis/2MASS/LGA/ 
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J n this Chapter we will examine the clustering dipole of galaxies from the Two Mi- 
cron All Sky Survey Extended Source Catalog as a function of increased depth of 
the sample. This analysis will allow us to constrain the /3 parameter from this dipole 
and consequently to estimate the density parameter Q^. Our method however is not 
to directly compare the peculiar velocity and acceleration of the Local Group; instead, 
we use the observed growth of the dipole to obtain these constraints. 

In the past, many different datasets have been applied to calculate the cluster- 
ing dipole. Generally speaking, there is no consistency on the amplitude, the scale of 
convergence of the dipole, and even on the convergence itself. The pioneering works 
used the revised Shapley-Ames catalog (Yahil et al., 1980) and the CfA cata- 
log (Davis & Huchra, 1982). A great advancement came with the launch of the 
far-infrared IRAS satellite and catalogs obtained thanks to this mission. The LG 
dipole from IRAS was studied first from two-dimensional data only (Yahil et al., 
1986; Meiksin & Davis, 1986; Harmon et al., 1987; Villumsen & Strauss, 
1987), then with redshifts included (Strauss et al., 1992; Schmoldt et al., 
1999; Rowan- Robinson et al., 2000; D'Mellow et al., 2004; Basilakos & 
Plionis, 2006) and with optical data added (Lahav et al., 1988; Kaiser & La- 
HAV, 1989). Samples with optical data only were also used (Lahav, 1987; Hudson, 
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1993), as well as galaxy clusters (Plionis & Valdarnini, 1991; Brunozzi et al., 
1995; KOCEVSKI & Ebeling, 2006). Among the most recent analyses one finds 
those directly related to the study presented here, which used the data from the Two 
Micron All Sky Survey. Maller et al. (2003) used the 2MASS Extended Source 
Catalog, concluded convergence of the clustering dipole from flux data only and used 
it to calculate the average mass-to-light ratio in the near-infrared Kg band and to 
estimate the linear biasing parameter hx-^ Erdogdu et al. (2006a) studied the 
acceleration of the LG from the 2MASS Redshift Survey, claimed convergence already 
at ~60 Mpc/Zi and estimated the (5 parameter by comparing the dipole with the LG 
velocity. In a more recent work, Lavaux et al. (2010) used an orbit-reconstruction 
algorithm to generate the peculiar velocity field for the 2MRS, extended it to larger 
radii, and observed no convergence of the clustering dipole up to at least 120 Mpc//i. 

Let us shortly reiterate some considerations of Section 1.5. Since both the gravita- 
tional force and the flux of a galaxy are inversely proportional to its distance squared, 
the clustering dipole can be calculated by assuming some average value of the mass-to- 
light ratio and summing only flux 'vectors' of galaxies (with positions on the sky used 
as angular coordinates). Such a dipole does not use distance nor redshift information 
of particular objects, which hinders the direct estimation of the depth up to which 
the dipole is measured and can bias the measurement of LG acceleration from it. This 
drawback is partially removed if a galaxy redshift survey (with redshifts of galaxies used 
as proxies of their distances) is used for the analysis; however, up to date the densest 
and deepest all-sky redshift surveys do not reach farther than to ~ 100 Mpc//i and 
contain no more than several dozen thousand sources. Owing to some recent claims 
of large-scale flows in the local Universe (e.g. Watkins et al. 2009; Kashlinsky 
ET AL. 2008), it is important to estimate the dipole from as deep catalogs as possible, 
even if they do not contain redshift measurements. We have thus decided to trade the 
advantages of redshift measurements for a huge number of galaxies and a much greater 
depth of a photometric-only catalog. The angular dipole thus obtained is additionally 



From now on we denote the general biasing parameter as h and the one in the Ks band as bx- 
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free from any redshift distortions, in particular from the rocket effect (Kaiser, 1987). 

In this Chapter, we focus on the issue of the convergence of the dipole, by analyzing 
its growth with increased depth of the sample. For that purpose we use the data from 
the 2MASS XSC, similarly to what has been presented in Maller et al. (2003). 
One of our goals is to expand and refine the latter work. We do it by modifying the 
criteria of removing Galactic objects from the sample (Section 3.3) and masking and 
filling the Galactic plane and bulge (the Zone of Avoidance, Section 3.4), but keeping 
in the ZoA the brightest and closest galaxies of big influence for the local motion (Sec- 
tion 3.5). We analyze the growth of the dipole and check how its direction changes on 
the sky. Next, we compare the observed growth of the dipole with theoretical expecta- 
tions, using the condition of known peculiar velocity of the Local Group, in a similar 
manner as was done in Lavaux et al. (2010). From this comparison we evaluate 
the P parameter and, by using externally provided value of bias, constrain the cosmic 
density Q^i- 

The Chapter is organized as follows. Section 4.1 focuses on the observed growth of 
the 2MASS clustering dipole. Next, Section 4.2 deals with the position of the 2MASS 
dipole on the sky. In subsequent Section 4.3 we discuss the results: first, in Subsection 
4.3.1 we present the theoretical basis to analyze the growth for known LG velocity; 
Subsection 4.3.2 focuses on the observational window of our measurement; finally, in 
Subsection 4.3.3 we compare our results with theoretical expectations and use this 
comparison to estimate the /3 parameter from the growth of the dipole alone. We 
shortly summarize and conclude this Chapter in Section 4.4, where we also provide the 
obtained value of the flm parameter. 

4.1 Observed growth of the 2MASS clustering dipole 

In this section we will use the data prepared as has been described in Chapter 3 to 
calculate the clustering dipole of the galaxies from the 2MASS XSC and analyze its 
growth. For that purpose we start with the formula for the flux dipole of the survey 
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given by Eq. (1.56) and change units of d into km/s, defining the scaled dipole as 

The flux of each galaxy is calculated from its magnitude nii as 

Si = So 10"°-^'"' , (±2) 

where 5*0 is the flux for a 0-magnitude object. As was already stated, we consider 
magnitudes in the Kg band, which was the main ('target') band of the 2MASS survey 
(for simplicity of notation, we sometimes skip the 's' subscript). Then Eq. (4.1) takes 
on the form: 

TT N 

d=fl^ S{0 mag) 10-°-^(^»+^PO)r, , {43} 
3k i 

where ^(Omag) = 1.122 x 10"^° Wm'^ (±1.685%) and the zero point offset ZPO = 
0.017 ± 0.005 (Cohen et al., 2003). The Ki magnitudes in Eq. (4.3) include also 
a negative offset of A = —0.2 mag added due to the underestimation of total fluxes 
by the isophotal magnitudes in the 2MASS XSC The quantity jx is the luminosity 
density in the K band. It is obtained for example from the integral (1.54) using 
the luminosity function in this band. The value of jx has been estimated by many 
authors in the recent decade: Kochanek et al. (2001); Cole et al. (2001); 
Bell et al. (2003); Eke et al. (2005); Jones et al. (2006); Smith et al. 
(2009); Lavaux & Hudson (2011) and the results from one work to another do 
not always agree within the errors quoted. Here, we use the determination that we 
consider the most reliable, calculated by Jones et al. (2006) from the luminosity 
function of more than 60,000 galaxies in the 6-degree Field Galaxy Survey (6dFGS, 
Jones et al. 2009): Jk = (5.9 ± 0.6) x lOHL^Mpc"^ (where ~ 6.8 x lO^^W, 
RiEKE ET AL. 2008^). Note that as the luminosity density depends linearly on the 
Hubble constant, the h factor cancels out and Eq. (4.3) may be rewritten as 

ax ^ 

d = C J2 lO-^ '^^^ri , {±4} 

Kynin 



^Lq is the Solar luminosity in the K band, calculated based on table 7 of Rieke et al. (2008), 
http: //www. iop.org/EJ/article/1538-3881/135/6/2245/aj 271287_mrt7.txt. 
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where C ~ 2620 km/s for parameter values as given above. The hmits of the sum- 
mation are now the minimum and maximum K magnitudes of galaxies in the sample. 
The lower limit in Eq. (4.4) is the magnitude of the brightest object in the catalog 
(excluding LG galaxies) and the upper one will be increased, as will be discussed later 
in the text. To retain the reliability of the analysis, we must note that the results 
remain trustworthy only up to the completeness limit of the catalog: the sample be- 
comes incomplete for objects with K > 13.6 mag (Jarrett, 2004). 

The growth of the dipole was calculated by incrementing i^max in the sum given in 
Eq. (4.4). Results are illustrated in Figures 4.1 and 4.2. Figure 4.1 shows the growth of 
the dipole as a function of the number of galaxies used for the calculation and presents 
also Galactic Cartesian components of the acceleration [d^ points towards the Galactic 
center, I = 0° , b = 0° ; dy is perpendicular to it in the Galactic plane; dz is perpendic- 
ular to the Galactic plane). The similarity of this plot to figure 1 of Maller et al. 
(2003) is intentional and allows us to draw the conclusion that although our data anal- 
ysis was slightly different to that of Maller et al. (2003) (altered shape and filling 
of the ZoA, different removal of Galactic objects), the two approaches are qualitatively 
the same. Nonetheless, we want to stress that such a presentation of the results may 
be misleading. First of all, as the discussed scales can be treated as Euclidean, for a 
uniform distribution the number of sources in a given radius is A^(r < i?) oc i?^, so a 
linear scale in the number of galaxies on the abscissa 'compresses' the left-hand side 
of the plot, while 'stretching' the right-hand one. This may suggest the convergence of 
the dipole, which was actually concluded by Maller et al. (2003). However, local 
inhomogeneities complicate this N{r < R) relation and finally a linear scaling in the 
number of galaxies up to a given flux threshold is not related in a simple manner to a 
scale expressed in distances nor in cutoff magnitudes. For that reason, what we need 
in order to draw proper conclusions about the growth and possible convergence of the 
dipole, is a linear scale in distance on the abscissa. 

Neither distances, nor even redshifts (including photometric ones) are currently 
measured for the whole 2MASS XSC (although some attempts are being made regard- 
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Figure 4.1: Growth of the 2MASS clustering dipole as a function of the number of galaxies 
used for the calculation (bottom axis), ordered by their Ks magnitudes (top axis). The 
thick red line is the amplitude of the dipole; the thin dotted and dashed lines (green, blue 
and magenta) are the Cartesian components (in Galactic coordinates). Two vertical lines 
illustrate the limits of the 2MRS Kg < 11.75 mag sample and completeness of the 2MASS 
XSC {Ks = 13.6 mag). Such data presentation could suggest that the 2MASS dipole has 
converged within sample limits, which is not the case. 



ing photo-2;'s, see Jarrett 2004 and Francis & Peacock 2010). We thus decided 
to define effective distances of galaxies based on their fluxes, with the use of the lumi- 
nosity function (LF) in the K band. If all the galaxies had the same luminosity, say 
L*, the relation between the observed flux S and distance r would be straightforward: 
r = ^L^/AnS. However, galaxies have different intrinsic luminosities, and obviously 
their LF is not a Dirac's delta (on the contrary, it is very broad). Therefore an esti- 
mated distance of a galaxy with a given flux must have a scatter. Constructing the 
estimator, a flrst choice could be the conditional mean, i.e. the expectation value for r 
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Figure 4.2: Growth of the 2MASS clustering dipole as a function of increased maximum 
Kg magnitude of the subsample (bottom axis). Corresponding effective distance is given at 
the top. The thick red line is the amplitude of the dipole; the thin dotted and dashed lines 
(green, blue and magenta) are the Cartesian components (in Galactic coordinates). Two 
vertical lines illustrate the limits of the 2MRS Kg < 11.75 mag sample and completeness of 
the 2MASS XSC {Kg = 13.6 mag). The lack of convergence of the dipole is clearly visible. 
Note also the steady growth of the amplitude as well as of the Galactic y component of the 
dipole. 

given S. Instead, we think that it is better to choose the conditional median for (a 
median value of distance given the flux). We consider it being more adequate to our 
problem: the same number of galaxies with a given flux have distances smaller and 
greater than the median, by definition. However, as we show in Appendix B, the two 
estimators differ by only a few percent. For the /C-band LF as given by JONES ET AL. 
(2006), i.e. a Schechter (1976) function 
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with = -23.83 + 5 log /i± 0.03 and a = -1.16±0.04, this effective distance for the 
magnitude K equals to 



(detailed calculations are presented in Appendix B; see also Peebles 1993). This 
proxy of distance is used in Figure 4.2, which differs from Figure 4.1 by a different 
scaling of the abscissa. The growth of the clustering dipole up to the completeness 
limit of the 2MASS XSC is now evident. Additionally, note that the growth has an 
essentially constant slope for r^fi > 150 Mpc/h, i.e. Kg > 12 mag. 

An interesting feature is the behavior of the Galactic Cartesian components of the 
dipole. The x and z ones are virtually constant for r^s > 150 Mpc/h; however, the 
y component still grows even at the limit of the catalog, similarly as does the total 
amplitude. This could point to some systematic effect, related to masking and filling 
of the Zone of Avoidance. We have however checked that the same qualitative behavior 
of the three components is observed for different shapes of the mask and the way it is 
filled; what is more, the effect exists even if we calculate the dipole having removed from 
the catalog all the galaxies with |6| < 10° (leaving the resulting strip completely devoid 
of galaxies). We thus conclude that the effect is real and must be due to the large-scale 
matter distribution in our cosmic neighborhood. An indirect confirmation comes from 
recent measurements of the local bulk flow, i.e. the volume-averaged peculiar velocity 
of galaxies in a several-dozen megaparsec sphere around us. Different authors (e.g. 
Watkins et al. 2009; Nusser & Davis 2011) who do not necessarily agree on 
the amplitude of the flow, still find consistent results when it comes to the direction 
(namely, Im — 280°, 6bf — 10°, i.e. roughly along the Galactic y-axis). The simplest 
explanation of the pull could be a single supermassive structure at a large distance, 
hidden behind the Galactic plane (as there is no known supercluster in this region just 
above the ZoA, see Figure 3.1). However, an alternative and perhaps more probable 
option is the summed effect of several sources, such as the Shapley Concentration and 
the Horologium Supercluster. This issue still awaits quantitative studies. 

We have also observed that adding to the sample the LGA galaxies that were not 



reff ~ 0.59 X 10°-^^ Mpc/h 




52 



4.2. MISALIGNMENT ANGLE 



present in our basic catalog has virtually no influence on the amplitude of the dipole 
and only slightly changes the values of particular components (see also below). 

4.2 Misalignment angle 

The linear theory relation (1.48) between the peculiar velocity and the acceleration 
of the Local Group predicts that the two vectors should be parallel. In reality, a non- 
zero misalignment angle between them is expected, due to various reasons (cf. Chap- 
ter 2). Observed values of this angle are usually of the order of 10° ~ 25° (Strauss 
ET AL., 1992; SCHMOLDT ET AL., 1999; Maller et al., 2003; Erdogdu et al., 
2006a) and we have obtained similar results, which confirms the validity of the linear 
approximation for the scales of interest. Figure 4.3 shows how the dipole direction on 
the sky changes as the sample depth increases from 'min', i.e. 5 Mpc/Zi, to 'max', equal 
to 310 Mpc/h (expressed here in terms of the effective distance). The black disk shows 
the CMB dipole. For comparison, we present also two other results of 2MASS data 
analysis: the 2MASS dipole with 'cloning' the ZoA, from Maller et al. (2003) 
(green triangle labeled M03) and the flux-weighted 2MRS dipole in the CMB frame of 
Erdogdu et al. (2006a) (magenta square, E06). We increase the Galactic longi- 
tude / from left to right for easier comparison with relevant figures in these two papers 
(respectively fig. 2 of M03 and fig. 7 of E06). The direction that we find at the limit 
of the sample lies somewhere in between those two results. 

In Figure 4.4 we plot the misalignment angle between the 2MASS clustering dipole 
and the CMB dipole direction, as a function of growing depth of the sample. The two 
central curves illustrate the effect of adding the LGA galaxies to the catalog, and in 
particular of the three bright galaxies located in the ZoA (Maffei 1 & 2 and Circinus). 
A specific 'tug-of-war' between these structures, located on almost opposite sides of the 
Galactic plane, results in a general raise of the angle by ~1°, which is well within Icr 
errors of the measurement (estimated via bootstrap resampling, see Subsec. 4.3.3 and 
Appendix C). As can be equally seen, the brightest (and presumably the closest, cf. 
Appendix A) galaxies have the most influence on the value of the angle. Note also that 
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Figure 4.3: Position of the 2MASS clustering dipole on the sky as a function of increased 
depth of the sample. The grid shows Galactic coordinates: longitude / (labels on bottom) 
and latitude b (left). The min/max labels refer to the minimum and maximum depth of the 
sample for which the dipole was calculated, equivalent to effective distances of 5 Mpc/h and 
310 Mpc/h, respectively. The black disk marked by 'CMB' is the direction of the peculiar 
velocity of the Local Group (/ = 272°, b = 28°). The green triangle labeled 'M03' shows 
the 2MASS clustering dipole of Maller et al. (2003): I = 263°, b = 40° (with cloning in 
the ZoA). The magenta square marked by 'E06' is the direction of the flux- weighted 2MRS 
dipole in the CMB frame of Erdogdu et al. (2006a): / = 245°, b = 39°. Blue labels 
indicate important structures in the local Universe: the Virgo Cluster {D ~ 17 Mpc), the 
Hydra Cluster {D ~ 47 Mpc) and the Leo Supercluster (z > 0.031). 
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Figure 4.4: Misalignment angle between the 2MASS clustering dipole and the peculiar veloc- 
ity of the Local Group as a function of increased maximum Kg magnitude of the subsample 
(bottom axis). The corresponding effective distance is given at the top. Two vertical lines 
illustrate the limits of the 2MRS Kg < 11-75 mag sample and completeness of the 2MASS 
XSC {Kg = 13.6 mag). The red solid and green dashed lines are shown to illustrate the effect 
of adding the Large Galaxy Atlas to the sample. Dotted lines show la errors (only for the 
XSC sample without the LGA, to avoid overcrowding of the plot). 



already for an effective distance as small as 50 Mpc/h the misalignment angle reaches 
a value of ~22° (~0.35 rad), with a minimum of 17° ± 5° for — 300 Mpc/h. 

One can propose several reasons for this misalignment. The first possibility are 
influences from extremely large scales. We cannot address this issue directly here; note 
only that the claims of Watkins et al. (2009) or Kashlinsky et al. (2008) 
of large-scale bulk flows with values much larger than expected in ACDM have been 
recently put in doubt by NussER et al. (2011), NussER & Davis (2011) and Os- 
borne ET AL. (2011): the matter is thus far from being settled. The second option 
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could be nearby unsuspected structures in the Zone of Avoidance. This issue was partly 
addressed in Chapter 2, where the influence of one such partially obscured object - 
the Local Void - was analyzed. As was shown, the improper filling of the ZoA in the 
LV region cannot bias the angle by more than 1°. Additionally, already in 2004 it 
was quite certain that thanks to surveys in various wavelengths (such as HI or X rays) 
there should be no unknown nearby large-scale structures hidden behind the Milky 
Way (e.g. Fairall & Lahav 2005). Finally, such a big value of the misalignment 
angle can be due to improper accounting of very local influences; we have tested this 
possibility by adding the LGA galaxies to our sample. 

Note however that our statistical assumption that the flux is a good proxy of the 
gravitational force, which may be appropriate for the whole sample, can fail for individ- 
ual galaxies. In particular, we have used a constant mass-to-light ratio for all galaxies, 
which is probably valid when averaged over many of them, especially since we use 
the Ks-ha,nd luminosity, which is known to be a better tracer of stellar mass than 
optical and other NIR wavebands. Still, the M/L^ ratio should vary somewhat with 
morphology and luminosity. For galaxies with smallest magnitudes, these variations 
may significantly affect their contribution to the direction of the clustering dipole. The 
same was noted by Maller et al. (2003), who also concluded that the 16° difference 
between their 2MASS dipole and the CMB one was caused by the brightest galaxies. 
We obtain a similar result here and interpret it as due to a specific non-linear effect: 
a few sources with the largest fluxes have a large influence on the misalignment. 

The high significance of the most luminous (in terms of the observed flux) galax- 
ies for the direction of the dipole confirms theoretical predictions of Chodorowski 
ET AL. (2008): the misalignment angle could be lowered if the observational window 
used was the optimal one, which would be the case if we removed those brightest galax- 
ies from the sample (as was also tried both by Maller et al. 2003 and Erdogdu 
ET AL. 2006a). This issue will be addressed in the following Chapter 5, concerning 
the optimized estimation of P from the 2MASS dipole. 
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4.3 Discussion and estimation of the p parameter 



Our results showing that the clustering dipole of 2MASS galaxies still grows even 
at the completeness limit of the sample, taken at face value, are consistent with what 
has been obtained by several authors, who used various datasets and methods: redshift 
survey of IRAS galaxies, Strauss et al. (1992); IRAS PSC^;, Schmoldt et al. 
(1999) and reanalysis by Basilakos & Plionis (2006); X-ray selected clusters, 
KOCEVSKI & Ebeling (2006); reconstructed velocity field of 2MRS, Lavaux et al. 
(2010). On the other hand, they contradict claims of convergence at scales even as 
small as 60 ~ 100 Mpc//i: optical sample of Hudson (1993); redshift sample of 
Abell/ACO clusters, Plionis & Valdarnini (1991); Brunozzi et al. (1995); 
IRAS PSCz, Rowan-Robinson et al. (2000); IRAS V^Gz and BTP, D'Mellow 
ET AL. (2004); 2MRS, Erdogdu et al. (2006a). Our analysis also suggests a 
different interpretation of the results of Maller et al. (2003) — data presentation 
as in our Figure 4.2 instead of figure 1 therein would possibly point to similar lack of 
convergence. On the other hand, in order to be able to directly compare our results 
with those of Erdogdu et al. (2006a), we would have to apply the same weighting 
as was done for the 2MRS sample, namely by the inverse of the flux-weighted selection 
function. We are unable to do it, not knowing distances nor redshifts for the whole 
sample. 

Apart from galaxy weighting, the discrepancies between the above listed results 
most probably stem from the different nature of catalogs and methods used for the 
calculation, and in particular may be due to distinct observational windows. Such 
a window for a given survey describes the sample: it may be interpreted as a filter 
(in real or Fourier space) through which we observe the Universe. Knowledge of the 
observational windows, necessary to correctly confront results as those given above, is 
also essential if we want to make comparisons with theoretical expectations. 
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4.3.1 Theoretical framework 

We would now like to check if the behavior of the 2MASS flux dipole is consistent 
with the predictions of the currently favored cosmological model, namely Lambda- 
Cold-Dark-Matter (ACDM) . We start by presenting the theoretical framework for such 
a comparison. It was first derived in the context of then-popular models like cold- 
dark-matter and isocurvature baryon, as described in detail in two classic papers: 
JUSZKIEWICZ ET AL. (1990) and Lahav et al. (1990). More recently, this approach 
was taken by Lavaux et al. (2010), who reconstructed the local peculiar velocity 
field (up to ~150 Mpc/h), applying the data from the 2MRS. The basic quantity for 
these comparisons is the joint probability distribution function for v and g, assumed 
to be a multivariate Gaussian (for more on this p.d.f., see next Chapter). 

In our case, we want to find the expectation value for the amplitude of the accel- 
eration of the Local Group knowing its peculiar velocity. As was already mentioned, 
the latter equals to i^^mb = 622 ± 35 km/s in the direction (/, b) = (272° ± 3°, 28° ± 5°). 
In the analysis of the present Chapter we do not use the directional constraint: the 
misalignment angle is integrated out (c.f. JuszKlEWicz et al. 1990). In the next 
Chapter, concerning the maximum-likelihood determination of /3 from 2MASS XSC, 
we will use the fact that the misalignment angle is known. Here, the theoretical pre- 
diction gives us the conditional velocity, (in units of km/s), which is related to 
the acceleration via the linear-theory Relation (1.55). The relevant formula for its 
amplitude has the form given by eq. (8a) of JuszKlEWicz et al. (1990) (note a typo 
therein) : 



V, 
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where we have used the following quantities: 

• Vm is the velocity induced by a single realization of the density field d{r) given 
the assumed power spectrum of density fluctuations: 
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with W{r) the window of the survey, included to mimic the velocity measured 
from a given survey; 

0",^ is the predicted (ensemble average) 1-D variance of the velocity measured 
from the survey: 



6 TT^ 



with P{k) = (|5fcP) being the aforementioned power spectrum, k the wavevector 
and w{k) the observational window in Fourier space; 

g is the cross-correlation coefficient of Vm and f : 



Gtt'^ a an 



Jp{k)w{k)dk, {aIo] 





where a is given by Eq. (4.9) with w{k) = 1. 

The second moment of the conditional velocity is given by eq. (8b) of JuszKlEWicz 
ET AL. (1990): 



hence the variance is equal to 



a2 = /i,-t;2. {4A2} 



The presence of the factor /(^m) in Eqs. (4.8)-(4.9) and the fact that we also have 
to include biasing b mean that in order to properly compare the conditional velocity 
Vc given by Eq. (4.7) with the amplitude of the dipole, d, measured from the data (Eq. 
4.1), the latter has to be rescaled, using a best-fit (3 parameter, as in Eq. (1.55). In 
fact, this rescaling enables us to estimate /3 in a straightforward way; we come back to 
this point a few paragraphs later. 
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4.3.2 Observational window 

In this subsection we discuss the proper observational window for our measurement. 
We start by noting that when calculating the peculiar gravitational acceleration of the 
Local Group from the dipole of the all-sky galaxy distribution, in general there are two 
schemes to postulate a relation between this distribution and that of the underlying 
mass (see e.g. Erdogdu et al. 2006a). The first scheme, called number weighting, 
assumes that the mass is distributed in the Universe as a continuous density field, 
which is sampled by galaxies in a Poissonian way. The second prescription, namely 
flux weighting, uses the assumption that all the mass in the Universe is locked to the 
mass of the halos of luminous galaxies. The 2MASS dipole as calculated in Section 4.1 
is of the flux- weighted type. 

The observational window W{r) of a fiux-limited survey is just its selection function 
ip{r). In the number-weighted scheme, this function measures what fraction of galaxies 
located in the distance interval (r, r + dr) are included in the survey: 

oo 

Vir) = , {±13} 



where ^{L) is the luminosity function and S^i^ is the limiting flux of the survey (e.g. 
Davis & Huchra 1982; Yahil et al. 1991). In our case of the fiux-weighted dipole, 
as given by Eq. (1.56), the appropriate window is the fiux-weighted selection function, 
\&(r), defined in such a way that 1 — \I^(r) is the percentage of light from a distance r 
which is not visible for the survey (Erdogdu et al., 2006a; Chodorowski et al., 
2008). It is given by 

oo 

/ L$(L)dL 

Ws{r) = ^-^ , {TU} 

/L$(L) dL 



where Lmin = 47rr^5'min. The window Ws{r), where the subscript S emphasizes the 
dependence of the window on the flux limit of the catalog, smoothly decreases to zero 
with increased distance from the observer. For a detailed discussion see Chodorowski 

ET AL. (2008). 
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On the other hand, 'ideal' surveys would be distance- (or volume-) limited: in this 
case there is no loss of signal up to a limiting distance R. Then the window function is 
simply unity for r < i? and zero otherwise (the so-called top-hat window). It has thus 
a form of a Heaviside step function; for a spherical survey we have 

WR{r) = Q^{R-r) [Hsj 

(e.g. JUSZKIEWICZ ET AL. 1990). In order to better reconstruct the dipole of the 
galaxy distribution from a flux-limited catalog, volume-weighting of the survey is com- 
monly mimicked by weighting individual galaxies by the inverse of the selection func- 
tion, at the expense of increasing shot noise from large scales (Strauss et al., 1992) 
and at a price of a possibly large rocket effect (Kaiser, 1987). The top hat is then the 
relevant window. In our case however, as we are dealing with angular (photometric) 
data only, we cannot estimate the selection function even in redshift space. Therefore, 
we do not weight galaxies in the sample (or in other words, we assign unit weights to 
them) and the relevant window function is given by (4.14). 

The Fourier-space counterpart of the observational window, w{k), is obtained for 
the clustering dipole from the real-space one as 

oo 

w{k) = kj W{r) ji{kr) dr (aI6] 



(Kaiser & Lahav, 1989; Lahav et al., 1990; Juszkiewicz et al., 1990), 
where ji is the spherical Bessel function of the first kind of order 1. Note that this 
window is not the Fourier transform of W{r). In our case of Ws{r) given by (4.14), 
the Fourier-space window reads (eq. 59 of Chodorowski et al. 2008) 

or, oo 

Wsik) = 1 - j^jool^l^^^ I sm(fcr) L^i, $(L^i,) dr , [417] 

where S^i^ and L^[^ are as above. The window Wg is a function of both the wavenum- 
ber k and the minimum flux Smin (hence the subscript s). On the other hand, it does 
not depend on the distance, as the latter is integrated out, together with absolute lu- 
minosities (with the use of the luminosity function). This means that when comparing 
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Figure 4.5: Comparison of tliree observational windows in Fourier space: the one of the 
2MASS flux-limited survey (thick red line) for a maximum magnitude of 13.6 mag and 
two for distance-limited surveys (top-hat) with iimax = 310 Mpc/h (blue thin line) and 
-Rmax = 133 Mpc/h (black dashed line). The former radius is the effective distance of 
galaxies at the limit of the 2MASS XSC; the latter corresponds to K = 11.75 mag, the limit 
of the 2MRS. 

observations to theoretical expectations, we never use the effective distance (4.6) for 
the calculations of the conditional velocity (4.7) and its variance (Tc (4.12). Our 
observable is the minimum flux of the sample, related to the maximum magnitude via 
'S'min = 'S'o 10"^-'^^"^^^ . For clarity however, in the relevant plots we prefer to show the 
results in terms of linear scaling in the effective distance. 

In Figure 4.5 we plot the 2MASS flux-limited window with /i'max = 13.6 mag 
and compare it with two top-hat windows for distance-limited surveys: one with 
-Rmax = 310 Mpc/h and the other with -Rmax = 133 Mpc/h, which are respectively the 
effective distances of galaxies at the limit of the 2MASS and 2MRS /'^max = 11-75 mag 
sample. The Fourier form of the top-hat window for the clustering dipole is (e.g. 
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JUSZKIEWICZ ET AL. 1990): 

with Jo = sinx/x being the spherical Bessel function of the first kind of order 0. The 
oscillating behavior of this window in Fourier space is the result of rapid truncation 
in real space. As we can see, the top-hat window with -Rmax = 310 Mpc//i, which 
asymptotically behaves in the same way as the 2MASS one with A'max = 13.6 mag, 
passes clearly more large-scale signal [k < 0.02 /i/Mpc) than the latter does. This 
will result in much faster convergence of the theoretical dipole measured through the 
top-hat window than through the 2MASS flux- weighted one (see below). On the other 
hand, the top-hat window with a cutoff equivalent to the limit of 2MRS at 133 Mpc//i 
blocks almost all the signal already for k < 0.01 /i/Mpc, as expected. 

4.3.3 Comparison with theoretical predictions 

We start by reminding that the analysis in this Chapter is performed within the lin- 
ear theory. For that reason, we do not include non-linear effects in the cross-correlation 
coefficient q (Eq. 4.10). Mathematically, this means that we set to unity two functions: 
the ratio of power spectra of density and scaled velocity divergence, TZik), and the co- 
herence function, C{k). In the next Chapter we relax these assumptions and the two 
functions will be discussed in greater detail there. 

Our goal here is to compare the observational data to the expectations of the 
ACDM concordance model. Hence the power spectrum used in Eqs. (4.9) and (4.10) 
is the spectrum of cold dark matter, with baryon effects included, as given by ElSEN- 
STEIN & Hu (1998). In this framework, P{k) oc k'^^T'^ik), where is the spectral 
index of primordial fluctuations and the transfer function T{k) depends on the param- 
eters h (Hubble constant divided by 100 km/s/Mpc), Vt^, f2b (density parameter of 
baryons) and Ug (present value of root-mean-square density contrast of matter fluctu- 
ations within a sphere of 8 Mpc//i). In our calculations we use the following set of 
parameters obtained from WMAP seven-year observations (LARSON ET AL., 2011): 
h = 0.71, Qrah"^ = 0.1335, n^^h^ = 0.02258, = 0.963 and as = 0.801. 
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Figure 4.6: Theoretically predicted growth of the conditional velocity of the Local Group for 
known v^^^ using two different observational windows: the one of 2MASS (blue dotted line) 
and top-hat (green dashed line), both calculated for the ACDM as the underlying cosmological 
model. The horizontal solid line is the observed velocity of the Local Group with respect to 
the CMB. For comparison, the prediction for the currently disfavored Standard CDM model 
with the top-hat window is shown with the black dot-dashed line. 



Before applying the observational data, in Figure 4.6 we compare the expected 
growth of the conditional velocity for the 2MASS window and for the top-hat case, 
both calculated with the above ACDM model parameters. As we can see, for such 
a flux-limited survey as 2MASS, the predicted velocity (blue dotted line) is far from 
converging to the limit of f^MB ~ 622km/s (horizontal line) even for r^s — 350 Mpc/h. 
In the top-hat case (green dashed line), the expected convergence is much faster, as 
could have been deduced already from the comparison of the windows, presented in 
Figure 4.5. However, it should be noted that even for all-sky catalogs that include 
redshifts, for which we can effectively model the dataset with a top-hat window (like 
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in the case of the 2MRS), the convergence of the dipole is not hkely before some 
200 Mpc//i, opposite to the results of Erdogdu et al. (2006a), where it is claimed 
that the contribution from structure beyond 6000 km/s(= 60 Mpc/Zi) is negligible. 
Figure 4.6 shows that for the latter distance, the conditional velocity for the top-hat 
window has reached less than 75% of its final value. On the other hand, for the 
2MASS window, Vc attains O.QSf^^^jg no sooner than for r ~ 470 Mpc//i (~14.5 mag 
in the Ks band, not shown in the plot), far beyond the completeness of the 2MASS 
XSC As an instructive case, in Figure 4.6 we additionally plot the prediction for the 
once-popular 'Standard CDM' model, currently strongly disfavored by observations 
{h = 0.5, = 1, l^b = 0.05, ris = 1, (Tg = 1), with the top-hat window (black 
dot-dashed line). Note how rapidly the dipole is expected to grow in this case. We 
suppose this may be one of the reasons why fast convergence of the observed dipole is 
still taken for granted by some authors also in the framework of the ACDM model. 

The most important result of this Chapter is presented in Figure 4.7. It shows 
the observed growth of the 2MASS dipole compared with the conditional velocity 
calculated from Eq. (4.7). As was already discussed at the end of Subsection 4.3.1, we 
rescale observational data by a best-fit parameter /3 = Q^^/b. The errorbars of the 
measured dipole were obtained from 1000-fold bootstrap resampling of the catalog (for 
details see Appendix C). The /3 parameter was fitted by a minimum-x^ procedure, 
including both observational errors (from bootstrapping) and the theoretical (cosmic) 
variance. Note that formally we are 'over-fitting', as the datapoints are correlated 
with each other, due to integral (cumulative) nature of the measured dipole. The 
analysis could in principle be redone for the differential dipole (shells in 'flux space'). 
For that purpose one would need to develop the theoretical part, which may not be 
straightforward especially for the second moment and so far has not been presented in 
any work that we are aware of. 

The observed growth, once rescaled, is well within la range of the theoretical 
prediction (apart from the datapoints for the smallest distances, where the number 
of galaxies is very small, hence the measurement is noisy). Remarkably, the best 
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Figure 4.7: Growth of the 2MASS clustering dipole (red points with la errorbars), compared 
with the theoretical expectation for the conditional LG velocity in the ACDM model (solid 
blue line with dotted blue lines showing la variances). The datapoints were rescaled by the /3 
parameter (see text). The horizontal solid black line is the Local Group velocity with respect 
to the CMB. Vertical lines give the Umits of the 2MRS and 2MASS XSC catalogs. 

agreement is obtained in the range 100 Mpc/h < < 300 Mpc/h, i.e. roughly 
between the extent of the 2MRS catalog and the 2MASS XSC completeness limit. The 
result of the fitting gives /3 = 0.38 with la confidence intervals of ±0.02. Owing to 
the considerations of Section 3.2 concerning the scatter introduced by the offset added 
when passing from isophotal fiducial elliptical aperture magnitudes to total fluxes, as 
well as due to findings of Chapter 2 on the bias introduced by the fact of possibly 
improper masking of the Local Void, we double this formal error. Hence our estimate 
is 

/3 = 0.38 ± 0.04 , [4I9) 
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precise within 11%. This measurement is in agreement with the analysis of Erdogdu 
ET AL. (2006a), where also data from 2MASS were used, although from much smaller 
scales. The clustering dipole was calculated there from the 2MASS Redshift Survey, 
then complete up to ~ 105 Mpc/Zi (maximum magnitude in the Kg band equal to 
11.25 mag). Their result, based on claimed convergence of the flux-weighted dipole, 
was /3 = 0.40 ± 0.09. We note that with our method, using only '2MRS-like' data, i.e. 
up to Ks = 11.25 mag, the best-fit value would be /3 = 0.35±0.06, which is smaller, but 
still consistent with our measurement from the whole sample and the one by Erdogdu 
ET AL. (2006a). a somewhat larger value was obtained by Pike & Hudson (2005), 
where the peculiar velocity field within 65 Mpc/h predicted from 2MASS photometry 
and public redshift data was compared to three independent peculiar velocity surveys 
based on Type la supernovae, surface brightness fluctuations in elliptical galaxies, and 
TuUy-Fisher distances to spiral galaxies. The best-fit from this comparison was /3 = 
0.49±0.04 [although close inspection of tables 1-3 from Pike & Hudson (2005) shows 
that these errors are possibly underestimated] . On the other hand, our value of /3 agrees 
within the errors with that obtained by Davis et al. (2011), who reconstructed the 
cosmological large scale flows in the nearby Universe using the SFI++ sample of Tully- 
Fisher measurements of galaxies and compared it with the whole sky distribution of 
galaxies in the 2MRS to derive 0.28 < (3 < 0.37 (68.3% confidence). Finally, in a recent 
paper, NUSSER ET AL. (2012) used the galaxy distribution in the 2MRS to solve for 
the peculiar velocity field, estimated absolute magnitudes of galaxies and constrained 
P by minimizing their spread. They found a value in excellent agreement with our 
determination: /3 = 0.35 ±0.1. 



4.4 Summary and conclusions 

In this Thesis, as the estimator of the LG acceleration we use the clustering dipole 
of the galaxies from the 2MASS Extended Source Catalog, which contains positions 
and fluxes of almost one million sources. As we want to reach as far from the Local 
Group as possible with the state-of-the-art observations, we have decided to trade the 
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advantages of redshift surveys for the huge number of objects and unprecedented sky 
coverage of the 2MASS XSC. This allowed us to measure the dipole up to 310 Mpc//i, 
which is more than 2 times farther than for the deepest existing all-sky survey of galax- 
ies with measured redshifts, the 2MASS Redshift Survey. The price to pay was our 
inability to weight galaxies with the inverse of their selection function (which would 
require redshifts to be known), so we could not directly mimic an 'ideal', distance- 
limited survey, with our data. On the other hand, the flux dipole does not suffer from 
redshift-space distortions, and in particular from Kaiser's rocket effect. 

We have focused mainly on the issue of convergence of the dipole, which has been a 
subject of debate for more than two decades now. We have shown that the flux dipole 
of the 2MASS XSC does not converge up to the completeness limit of the sample, which 
is ~ 13.6 mag in the Ks band. Moreover, we find that beyond an effective distance 
of some 150 Mpc//i, the growth is induced mainly in the negative Galactic Cartesian 
y direction. We have checked that this behavior is not due to systematic effects re- 
lated to masking out and artificially filling the Zone of Avoidance. We thus suppose 
that the effect is related to large-scale matter distribution in our cosmic neighborhood. 
Additionally, the misalignment angle between the measured dipole and the velocity of 
the LG is found to be of the order of 20°. We also find that the misalignment comes 
mainly from some nearby galaxies, bright in the near infrared. This is a hint that 
their removal could optimize the analysis, as was already theoretically investigated by 
Chodorowski et al. (2008). We come back to this issue in the next Chapter. 

The observed lack of convergence of the dipole reflects the fact that for such a sur- 
vey as the 2MASS XSC, we measure the growth using only angular and photometric 
data. As the luminosity function of galaxies is very broad, the information coming 
from a given distance is 'spread' onto a wide interval of measured fluxes. However, by 
comparing this observed growth with theoretical expectations for known v^^^^, in the 
framework introduced by JUSZKIEWICZ ET AL. (1990) and Lahav et AL. (1990), 
we find consistency with the predictions of the ACDM model, once the proper obser- 
vational window of the 2MASS flux-limited catalog has been included. On the other 
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hand, in the case of the top-hat window, which is apphcable for a perfect, distance- 
hmited survey, the theoretically expected convergence of the dipole is much faster, 
and takes place at a scale of ~ 200 Mpc/h. This is considerably larger than the 
60 Mpc/h claimed by Erdogdu et al. (2006a), where galaxies from the 2MASS 
Redshift Survey were weighted by the inverse of their selection function to mimic such 
a volume-limited catalog. 

For the 2MASS XSC, we do not have redshifts for the whole sample and cannot 
define the selection function. Our inability to weight galaxies and to mimic a distance- 
limited survey, together with the lack of convergence of the dipole, do not allow us to 
measure the parameter /3 from a direct v^^ - g^^ comparison. However, we were still 
able to estimate this parameter due to the consistency of the observed growth with 
theoretical predictions; the result is /3 = 0.38 ± 0.04. This value is in agreement with 
and has better precision than the one from the analysis of Erdogdu et al. (2006a), 
who found (3 = 0.40 ± 0.09 from the clustering dipole of the 2MRS, which has a depth 
of ~ 100 Mpc/h. 

Having estimated the value of /3, we can go further and derive the density parameter 
D.jji- By definition of /3, we have f2m = (/3 6)^''°'^^, we can thus express our measurement 
as 

= 61-82 X (0.17 ±0.03) . {4~20] 

We can also constrain alone, provided that we know the linear bias, b. The latter is 
dependent on the wavelength and was estimated for the Kg band by Maller et al. 
(2005), who calculated the angular correlation function of galaxies in 2MASS and 
inverted it using singular value decomposition to measure the three-dimensional power 
spectrum. A fit of CDM-type power spectra in the linear regime allowed them to 
constrain the linear bias as = 1.1 ± 0.2. Using this value, we obtain the following 
estimate: 

= 0.20 ±0.08 , 

with a relative error of ~40%, much larger than for /3 alone due to the big uncertainty 
in the biasing parameter. Still, this result, supporting the fact that the density of non- 
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relativistic matter in the Universe is well bellow the critical one, roughly agrees with 
other independent dynamical estimates of Q^a, such as e.g. that by Feldman et al. 
(2003), who used mean relative peculiar velocity measurements for pairs of galaxies 
and obtained fim = 0.30lo;o7. Our findings are also perfectly consistent with the results 
derived from galaxy cluster data. For instance, Mantz et al. (2010) used obser- 
vations of the growth of massive. X-ray flux-selected galaxy clusters and derived, for 
spatially fiat models with a constant dark energy equation of state, = 0.23 ± 0.04. 
Last but not least, the density parameter derived here agrees within the errors with the 
one estimated from the analyses of the angular power spectrum of the cosmic microwave 
background temperature anisotropies, for instance those from WMAP seven-year ob- 
servations (Larson et al., 2011): = 0.267 ±0.029. 

The conclusions of this Chapter can be verified if new data become available. In 
particular, one could examine the growth of the 2MASS dipole with redshifts of galax- 
ies as their distance estimates, using proper weights. In the near future we cannot 
however hope for a uniform and sufficiently deep sample of spectroscopic redshifts for 
the catalog, even if the 2MASS Redshift Survey is continued (as it was planned to 
reach Kg = 12.25 mag, see Huchra et al. 2005). The calculation may be thus 
feasible only for photo- 2;'s, if they are available for the whole 2MASS XSC Some first 
steps towards this goal have been already taken (Jarrett, 2004) and a promising 
direction may be to cross-correlate 2MASS with the data from the Wide-field Infrared 
Survey Explorer (WISE, Wright et al. 2010) that are already available (Cutri 
ET AL., 2012). We can also take a different approach to estimate /3 from LG velocity- 
acceleration comparison. By maximizing the correlation between observed velocity and 
acceleration, we can use the maximum likelihood method to try to obtain an optimal 
measurement of /3. Such study is presented in the next Chapter, where we use the 
same 2MASS data but refine the method of statistical analysis. 
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Optimized estimation of parameter 



he analysis presented in Chapter 4 shows that the growth of the 2MASS cluster- 



ing dipole is consistent with ACDM predictions. This consistency allowed us to 
constrain the (5 parameter by comparing the observed growth with the theoretically 
predicted one. However, it is still possible to optimize our measurement. First of all, 
owing to the integral nature of the dipole, the points in Figure 4.7 are strongly corre- 
lated. Fitting a theoretical curve gives then a minimum much less than 1. Secondly, 
the window function Ws{r) as given by Eq. (4.14) is not the optimal one, which results 
in a significant value of the misalignment angle between the peculiar velocity of the 
LG and the 2MASS clustering dipole. As was shown e.g. by Strauss et al. (1992) 
and followed by Chodorowski et al. (2008; hereafter 'C08'), in order to obtain 
the best estimates of cosmological parameters by comparing the LG velocity and ac- 
celeration, one should preferably use the maximum likelihood approach and optimize 
the window through which the clustering dipole is measured. As one of the indications 
of this optimization, the misalignment angle is then expected to decrease. 

In this Chapter we will confront the theoretical predictions of COS with observa- 
tional data, optimize the window and use the maximum likelihood method to constrain 
the parameter (5. The Chapter is organized as follows. Section 5.1 presents the an- 
alytical description of the likelihood function for parameter estimation from the LG 
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velocity-acceleration comparison. In Section 5.2 we apply 2MASS data to the theoreti- 
cal framework of COS in order to determine the optimal window for the measurement. 
Section 5.3 deals with the optimized measurement of the /3 parameter. Finally, in 
Section 5.4 we summarize the Chapter. 

5.1 Analytical description of the likelihood 

Let us come back to the linear-theory Equation (1.48) relating the peculiar velocity 
of the LG with its acceleration: 

For a number of reasons, such as nonlinear effects, observational windows through 
which the velocity and acceleration of the LG are observed or shot noise from the 
edges of the catalog, this linear relationship is spoiled. Therefore, even if we had 
a perfect estimator of g^^, in order to measure /3 from this relation we could not 
simply equate the two dipoles, but a more sophisticated approach would be needed. 
Commonly adopted is the maximum likelihood approach, which provides a way to 
account for these effects and to compute errors of estimated cosmological parameters. 
We will now shortly sketch the analytical description of the method. For details see 
C1ECIEL4G & Chodorowski (2004), denoted hereunder as 'CC04', or COS. Note 
that for simplicity we will express both v^^ and g^^G units of km/s and as the 
theoretical quantity for the former (and its divergence) we will use the scaled velocity 
v^c = f~^{^m)Vohs, where / is the growth factor. We will come back to the 'observed' 
(measured) velocity of the LG when dealing with the expression for the likelihood. 
Finally, we will often refer to the LG acceleration as 'LG gravity' (or simply, 'gravity') 
and skip the subscript 'LG' in both gravity and velocity. 

Let p{g,v) denote the joint probability distribution function (p.d.f.) for the LG 
gravity and peculiar velocity. In a Bayesian approach, we ascribe a priori equal prob- 
abilities to values of unknown parameters, which allows us to express their likelihood 
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function, given v and g of the LG, via the p.d.f. for v and g: 

£(param.|g, v) = p{g, D|param.) . 

It is a standard practice to approximate the p.d.f. by a multivariate Gaussian 
(Strauss et al., 1992; Schmoldt et al., 1999) and this assumption was al- 
ready made implicitly in Section 4.3. Numerical simulations (Kofman et al., 1994; 
C1ECIEL4G ET AL., 2003) show that non-gaussianity of fully nonlinear g and v is in- 
deed small. This is not surprising since, e.g. gravity is an integral of density over a large 
volume (Eq. 1.41), so the central limit theorem should at least partly be applicable. 
Using statistical isotropy of g and v, their joint p.d.f. can be simplified to the form: 

x'^ + y^ — 2gfixy 



(2vr)3a3,g 



(53} 



2(1 - g^) 

(JUSZKIEWICZ ET AL., 1990; Lahav ET AL., 1990), where (Tg and are the theoret- 
ical r.m.s. values of a single Cartesian component of gravity and velocity, respectively. 
From isotropy, a"^ = {g ■ g)/3 and al = {v ■ v)/3, where (■) denote the ensemble av- 
eraging. Next, {x,y) = {g/crg,v/aj,), and fi = cos6, with 6 being the misalignment 
angle between g and v. Finally, g is the cross-correlation coefficient of Qm with v^, 
where Qm and Vm denote an arbitrary Cartesian component of respectively g and v. 
This coefficient was already introduced for a specific case in Eq. (4.10). Here we will 
provide its general form. From isotropy, 

{g ■ v) , , 

Also from isotropy, 

where Smn denotes the Kronecker delta. In other words, there are no cross-correlations 
between different spatial components. 

For a given all-sky galaxy survey, the LG gravity is measured effectively through 
the window of the survey, Wg (cf. Subsection 4.3.2): 

g = J^Sir)Wgir)r-^. {51} 
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In contrast, the LG velocity is not estimated from a velocity survey (i.e., from a catalog 
of peculiar velocities of galaxies), but measured directly from the dipole anisotropy of 
the CMB. Still, to relate it to theoretical quantities, we write: 



The quantity ^9 = — V ■ v is the minus scaled velocity divergence and we assume 
that the velocity field is irrotational. Kelvin's circulation theorem assures that the 
cosmic velocity field is vorticity-free as long as there are no shell crossings. N-hody 
simulations (e.g. PiCHON & Bernardeau 1999; Pueblas & Scoccimarro 2009) 
show that the vorticity of velocity is small in comparison to its divergence even in the 
fully nonlinear regime. Thus, similarly to also v can be expressed as a Coulomb 
(Newton) integral over its source, i.e. the field of the velocity divergence. Here we 
do not assume that we know the latter from observations, but we know from theory 
its statistical relation to the density field. This is sufficient for our purposes in this 
Chapter. Since v is directly measured from the CMB dipole, the effective velocity 
window, W^, which we have introduced in Equation (5.7), is essentially unity. This 
was implicitly assumed in Section 4.3.1. Here we modify slightly this form of the 
window to reflect the finite size of the LG. We adopt 



which has a small-scale cutoff, = 1 Mpc//i (van den Bergh, 2007; van der 
Marel & GUHATHAKURTA, 2008). A specific form of the gravity window, Wg, of 
the 2MASS survey was given in Subsec. 4.3.2. We will come back to the issue of its 
optimization later. 

In Fourier space, relations (5.6) and (5.7) read: 






1, otherwise , 



r < r, 




ik 



SkWg{k) 




ik 
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where the subscript k denotes the Fourier transform and w are windows in Fourier 
space. From relation (4.16) we obtain for the velocity window 

where j'o represents the spherical Bessel function of the first kind of order 0. 
From equations (5.9) and (5.10) we have 

and 

1 roo 

{vv) = — wl{k)n{k)P{k)dk , (513} 
Ziv^ Jo 



with 



being the ratio of the power spectra of scaled velocity divergence and of density. In 
Chapter 4 we have assumed this ratio to be unity; here, as we want to include possible 
non-linear effects, we relax this assumption. In particular, we will use the results 
obtained by CC04 based on numerical simulations. They found the following fit: 



7^(fc) = [1 + (7.071A;)^]-", [ 5. 15 J 

where 



a 



-0.06574 + 0.29195 ag for 0.3 < ag < 1 , [sle] 



with (Tg being the present value of the root-mean-square density contrast of matter 
fluctuations within a sphere of 8 Mpc/h. CC04 argued that the ratio of the power 
spectra practically does not depend on the background cosmological model. This ratio 
is unity in the linear regime {k <^ 1) but decreases in the nonlinear regime, because 
the velocity grows slower than it would be expected from the linear approximation. 
Furthermore, 

= ^ Wg{k) w^{k) C{k) Pl'\k) P'/\k) dk, [5l7] 
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where C{k) is the coherence function (CF), or the correlation coefficient of the Fourier 
components of the gravity and velocity fields: 

^ ' {\gk\'y/'{\vk\Y/' ilSklY^mW^' ^—^ 

(Strauss et al., 1992). Again, as in the case of the ratio of the power spectra, 
deviations of this function from unity were neglected in Section 4.3. It is however 
not constant and, similarly as was in the case of TZ, depends on the wavevector k 
and the erg parameter. This function was analyzed in detail by Chodorowski & 
C1ECIEL4G (2002) and refined by CC04, where the following fit was found, also 
based on numerical simulations: 

C{k) = [1 + {aok - a^k^-^ + aik^^] , {sJ9] 

with the coefficients given by scaling relations in erg: 

ao = 4.908(7^-^^°, 
ai = 2.663 (T^■^=^^ 
02 = 5.889 cr^-^^^ 

The fit was calculated for k G (0,1) h/Mpc and erg G (0.1,1), with the imposed 
constraint C{k = 0) = 1. This constraint assures that for sufficiently large, linear 
scales, the relation between the gravity and the velocity is deterministic and linear 
(Eq. 5.1). Chodorowski & Cieciel4G (2002) investigated numerically also the 
dependence of the CF on and found it to be extremely weak. 

Hence, we obtain finally for the cross-correlation coefficient of velocity and gravity: 

j;^wg{k)w^{k)c{k)n^/\k)p{k)dk 



1/2 



j,^wiik)P{k)dk [!^wi{k)n{k)P{k)Ak] 



1/2 



[5.21 



Equations (5.12), (5.13) and (5.21) specify all the parameters (the variances and the 
cross-correlation coefficient) that determine the joint p.d.f. for g and Equation (5.3), 
in the absence of observational errors. The deviation of the cross-correlation coefficient 
from unity is then due to different windows through which the gravity and the velocity 
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of the LG are measured, and due to nonlinear effects. The latter are described by the 
two functions: the CF and the ratio of the power spectra. Later on we will describe 
the corrections to the cross-correlation coefficient if observational errors are non-zero, 
as is in the case of our analysis. As in the preceding Chapter, when calculating the 
theoretical quantities that require the power spectrum of density perturbations, we use 
the framework of ElSENSTEIN & Hu (1998) with cosmological parameters derived 
from WMAP seven-year observations (Larson et al., 2011): h = 0.71, fim/i^ = 
0.1335, n^h^ = 0.02258, = 0.963 and as = 0.801. 

The final quantity necessary to calculate the joint p.d.f. of LG velocity and gravity 
is the observational window of the latter. It was derived and discussed in detail by 
C08. In Chapter 4 we have already presented its form when all the galaxies brighter 
than some limiting flux S'mm are included. However, as was shown by C08, in order to 
optimize the measurements which use the p.d.f., one should remove from the dataset 
the objects brighter than an appropriate minimum magnitude -ft'min- The reason for 
that procedure is that we need to mitigate nonlinear effects as much as possible. For 
large Fourier-space wavenumber k, the coherence function of velocity with gravity 
(Eq. 5.19) drops significantly below unity, decreasing the value of the cross-correlation 
coefficient (Eq. 5.21). Suppressing then the gravity window for large k has almost no 
effect on the cross-term (which is the numerator of Eq. 5.21), while it decreases the 
gravity variance, the square root of which appears in the denominator of this equation. 
This manipulation on the gravity window helps therefore to achieve the best possible 
correlation between the LG velocity and gravity. However, when one suppresses the 
gravity window for scales which are linear enough so that the CF is close to unity, 
one worsens the correlation again. This is so because even for linear fields (CF and 
the ratio of power spectra equal to unity) the cross-correlation coefficient decreases 
for increasingly different windows of velocity and gravity. As a result, for some value 
of Kmin (or equivalently, a maximum flux Sma.x oc lO"^-^^"""), at which we should cut 
the sample by removing bright galaxies, the cross-correlation coefficient will have a 
maximum. 
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The observational window of the 2MASS survey including all the galaxies between 
the maximum flux Smax and the minimum one Smm (the latter given now by the 
completeness of the catalog) is 

-^max 

/ L$(L)dL 
J L<^{L) dL 



(eq. 47 of COS) where Lmax = 47rr^S'max and ^{L) is the luminosity function of galaxies. 
In Fourier space this translates into 



I smikr)L^,^ $(L„i,) dr 



k!^L^{L) dL^ 

(eq. 54 of COS). Note however that the flux (or gravity) window does not appear 
in Equation (5.13) for the velocity variance, and in Equation (5.21) for the cross- 
correlation coefficient it appears in such a way that its absolute normalization cancels 
out. In other words, the p.d.f. is sensitive only to the shape of the gravity window. 

COS calculated the cross-correlation coefficient as a function of the minimum mag- 
nitude i^min for the window (5.23) and with the LP taken from Bell et al. (2003). 
They found that q has a maximum (~ 0.97) for i^min = 4.5 mag, under the assump- 
tion of negligible shot noise (from dilute sampling by distant galaxies of the underlying 
mass density field). Additionally, if the shot noise can be neglected, the minimum of 
the cross-correlation coefficient was shown to exactly correspond to the minimal ex- 
pectation value of the misalignment angle between the LG velocity and acceleration. 
On the other hand, should the shot noise prove to be significant, the optimal value of 
Kmin was predicted to increase. 

In order to obtain the final expression for the likelihood of cosmological parameters 
as a function of observables and theoretical quantities, based on Relation (5.2), we must 
account for observational errors. Let us then assume that for the LG acceleration they 
are given by a 3D vector, e. The variance of a single spatial component of measured 
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gravity, a^, is a sum of the ID cosmo logical component, a^^^, and errors, Because 
the gravity is inferred from the observations of a galaxian, rather than mass, density 
field, we must include biasing: 0"^^ = b'^{g^)/3. Thus, 



9,c 



As for the measured velocity of the LG, known from the CMB dipole, we follow 
CC04 and COS and neglect its observational errors for the purpose of our analysis. 
Its ID theoretical (cosmic) variance can be then calculated by plugging Eq. (5.15) 
into the integral (5.13). There is no biasing in this case, but in order to compare it 
with observables, we return to 'physical' variables, inserting the growth factor /(f2m) 
squared. The ID velocity variance is thus a function of VL„^\ 

where (f^)^^^ = 988 km/s is constant (depends only on the velocity window and the 
power spectrum, which are fixed). 

As was shown by CC04, including the observational errors has the effect of lowering 
the value of the cross-correlation coefficient. The effective cross-correlation coefficient 
will have the form 

(CC04; COS), where g is given by Equation (5.21). 

We can now proceed to the expression for the likelihood function. The observables 
('data') are the measured values of the LG velocity and gravity, and g^, respectively, 
as well as /im, i-e. the cosine of the misalignment angle (where g^ and /im vary when 
the window is changed). The model parameters ag and g depend on the biasing b, 
the gravity errors e and the gravity window Wg] a„ depends only on fim. Because of 
the relation ^l^^ = Pb, the minus log-likelihood function can be cast to the following 
form, depending explicitly on the biasing and /3: 



ln£ = 31n/3 + 31n[6aga„(l-^2)^/2j ^ 

1 f gl, . '^gjJ-mgmVm' 



2(1-^2) ^^2 ^252^2 f3baga^ 



const 
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(CC04; COS), where cr^ = (v^) /3 is the theoretical ID variance of scaled velocity. In 
order to obtain the most likely values for particular parameters, one has to find the 
minimum of the function given above with respect to them, given the observational 
data in hand. 

5.2 In quest for the optimal window 

We will now apply the 2MASS data and try to find the optimal window for our 
measurements. The catalog was prepared and analyzed in the same way as described 
in Chapter 3. The difference with respect to the analysis of the growth of the flux 
dipole presented in Chapter 4 is that now the measured dipole has a constant cutoff at 
the high- magnitude side, -ft'max = 13.6 mag, reflecting the completeness of the sample, 
and it is the minimum magnitude of the sample that is increased: 

where C ^ 2620 km/s as explained in Section 4. The amplitude of this dipole for 
a given f^min should be understood as the 'measured gravity' referred to in this 
Chapter. Note that removing the brightest (in terms of received flux) galaxies from 
the sample means that the estimated amplitude of the dipole will usually be smaller 
than the one measured from the whole dataset. However, this loss of signal will be 
compensated for by effectively 'weighting' gravity by its theoretical r.m.s. (5.12) with 
a relevant window. 

We will now try to find the optimal value of -ft'min for which we can make the 
maximum-likelihood estimation of (3. We start by examining the theoretically pre- 
dicted cross-correlation coefficient, g, as given by Eq. (5.21), which would be applicable 
in the absence of observational errors. It depends on the shape of the gravity window 
and we plot it as a function of increased -ft'mm in Figure 5.1. The cross-correlation 
coefficient reaches its maximum fPmax = 0.971 for -fTmin = 5 mag (at an effective dis- 
tance of 5.9 Mpc/h). This magnitude is somewhat larger then the i^min = 4.5 mag 
predicted by COS and in our opinion the main reason for that slight discrepancy is the 
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Figure 5.1: Comparison of the theoretically predicted cross-correlation coefficient g between 
LG gravity and velocity (blue dashed line, right-hand-side scaling) with the observed mis- 
alignment angle between LG velocity and the 2MASS clustering dipole (red points with 
errorbars, left-hand-side scaling), both as a function of increased minimum magnitude of the 
sample, and for a constant maximum one, -f^max = 13.6 mag. The top axis gives effective 
cut-off distances related to -ft'min- 



form of the luminosity function used to calculate the gravity window (5.22). COS used 
the Bell et al. (2003) LF, i.e. a Schechter (1976) function with a = -0.77 and 
M* = -23.29 5 log /i, while we take the updated LF by Jones et al. (2006), where 
a = —1.16 and M* = —23.83 -|- 51og/i. Note however that our value of Qma^ itself is 
very close to the one found by COS. 

COS emphasize that statistically, the maximum of the cross-correlation coefficient 
corresponds to the minimum of the misalignment angle (m.a.) provided that observa- 
tional errors, such as the shot noise, are negligible. However, owing to a broad p.d.f. of 
the expected angle (cf. figure 3 of COS), the scatter around this relationship is quite 
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significant. As can be seen in Figure 5.1, indeed if we vary tlie observational window 
(increase tlie minimum magnitude of the sample), the m.a. gradually decreases from 
the original ~ 19° to a minimum value. However, the smallest m.a. is attained for i^mm 
significantly larger than 5 mag. Analysis of the data shows that we have a relatively 
wide 'plateau' of small values of the angle with 6'min = 9° ±2° for -ft'min around 6.3 mag 
(rcff = 10.8 Mpc//i). Applying such a cutoff would be equivalent to removing some 50 
brightest galaxies from the sample. The la errorbars here, as well as in Figure 5.1 and 
hereafter, were estimated from 1000-fold bootstrap resampling of the catalog, similarly 
to what was done in Chapter 4 (see Appendix C for details). 

The cross-correlation coefficient for -ft'mm = 6.3 mag, where we have the minimum 
of the misalignment angle, 6'min = 8.9°, is p|6»min = 0.964. It is interesting to check how 
likely it would be in ACDM to observe such a misalignment angle, or smaller, for this 
value of the cross- correlation coefficient. For that purpose, we follow COS and use the 
conditional distribution of the m.a. with the LG velocity as a constraint (JuszKlEWicz 
ET AL., 1990; Lahav et al., 1990). In the small-angle approximation, the condi- 
tional m.a. is Rayleigh-distributed: 



e ( e 



2 



p(^|^)^-exp -— , {b^ 
with the characteristic measure of misalignment 

QV 

where y = v^^/a^. Integrating the distribution (5.29) from 6^ = 0° to 6'min = 8.9° with 
fixed p = 0.964, we obtain that observing such a misalignment angle, or smaller, at 
this magnitude cut has more than 56% probability. For Kmin = 6.3 mag, the observed 
m.a. is almost the same as its expectation value: 



\v) = {531} 



(COS), which in this case is equal to 8.74°. 

Coming back to the expression for the log-likelihood, Eq. (5.27), in order to calculate 
P, we will need to apply the effective cross-correlation coefficient, g. It includes gravity 
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Figure 5.2: Comparison of tlie effective cross-correlation coefficient between LG gravity and 
velocity, g (blue dashed line for bx = 1-1 and two dotted lines for bx = 0.9 and bx = 1-3, 
top and bottom respectively) with the cosine of the observed misalignment angle between 
LG velocity and the 2MASS clustering dipole (red points with errorbars), both as a function 
of increased minimum magnitude of the sample, and for a constant maximum one, ii'max = 
13.6 mag. The top axis gives effective cut-off distances related to i^min- 



errors (Eq. 5.26), but depends also on the r.m.s. of gravity (hence its observational 
window), as well as on the biasing parameter, (although very weakly). These facts are 
illustrated in Figure 5.2. The blue curve shows ^ for 6/^ = 1.1. The lower and upper 
dotted ones are respectively for bx = 0.9 and bx = 1-3. These values span the la 
interval of bx as determined by Maller et al. (2005). As can be seen, the value of 
biasing has a very small influence both on the position and on the value of ^max- What 
is more, there is a wide range of -ft'min for which the effective cross-correlation coefficient 
is virtually constant (for a given bx) and maximal. For example, if we take bx = 1-1, 
we have gma.x = 0.967 for -ft'min ^ (4.8,5.2). Situation is similar for other reasonable 
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values of h^- The effective cross-correlation coefficient is slightly lowered with respect 
to the theoretical one, although its maximum is not shifted towards bigger i^min- In 
the plot we also show the cosine of the measured misalignment angle, /im = cos6'ni, as 
this is the quantity that enters into the expression for the likelihood, Eq. (5.27). We 
can see that as was in the case of the theoretical cross-correlation coefficient, the -R'mm; 
for which the effective coefficient attains its maximum, is also below 6.3 mag, where 
/im is the biggest. We cannot thus unambiguously select a single window that would 
be optimal, based on this comparison alone. However, already at this stage we can 
support the claims of COS that -R'min = 8 mag is too big a magnitude at which one 
should limit the sample from below to optimize the window. Such a cutoff was tried by 
Maller et al. (2003), who claimed that removing the 375 brightest galaxies from 
their sample with < 8 mag resulted in a substantial decrease in the misalignment, to 
about 5.2°. This was then tested by Erdogdu et al. (2006a) from 2MRS data with 
no confirmation. We cannot confirm this result from our data either. If we remove 
the sources brighter than i^' = 8 mag (453 galaxies in our case^), the misalignment 
angle is equal to ~ 11°, which is indeed less than the ~ 19° obtained from the whole 
sample, but still significantly larger than 5.2°. What is even more important for the 
likelihood analysis, both the 'pure' cross-correlation coefficient and the effective one 
are for /i'min = 8 mag several per cent below their maximum values. These two effects 
will translate into the instable behavior of /3 estimated for such minimum-magnitude 
cuts. We present it in the following Section. 



5.3 Precise measurement of the p parameter 

The analysis presented in the previous Section did not provide a single value of 
Krain for which the window (5.22) may be unambiguously claimed optimal. We can 
however propose -ft'min — 6 mag as a 'reasonable' value here (the misalignment angle is 
close to its minimum and the effective cross-correlation coefficient is almost maximal). 

^This number is slightly different from the 375 given by Maller et al. (2003), but note that it 
includes LGA galaxies and also the 'clones' in the ZoA. 
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Such a minimum-magnitude cut would mean removing the 35 brightest galaxies from 
the sample and is equivalent to an effective distance of 9.4 Mpc//i. For the interested 
reader, the properties of these galaxies, such as their near-infrared magnitudes, mor- 
phological types, distance estimates and radial velocities, can be found in Appendix 
A. As can be seen in Table A.l, most of them are located at a distance smaller than 
~11 Mpc (~ 8 Mpc//i), so indeed removing them from the sample is roughly equiva- 
lent to including in the analysis only the sufficiently distant ones, thus mitigating large 
non-linear influences from those nearby. 

We can also calculate the (3 parameter for different cutoffs of the window. We 
should expect (3 to stabilize for an interval of minimum magnitudes, in a similar way 
what was observed for both q and /im- Let us then come back to the expression for the 
logarithmic likelihood, Eq. (5.27). We calculate its partial derivative with respect to /3 
and equate it to zero, to flnd the maximum of ln£. This gives the following quadratic 
equation (CC04; COS): 

b ag ¥ 

The relevant solution for /3 has the form: 

^MLE = ^^TT — (\/^^/^m^m + 12(1 - ^2)^2 _ ^^^^^ ^ [533] 

where the effective cross-correlation coefficient, is given by Eq. (5.26) and the ID 
variance of gravity dg is calculated via Eq. (5.24). 

The important thing to note is that the dependence of Pmle on h is more compli- 
cated (and much weaker) than in the case of the 'simple' /3 cx relation, because both 
Q and dg are functions of h. This is readily visible in Figure 5.3, which shows how /3mle 
chang function of increased minimum magnitude of the sample, for three values 

of the /C-band biasing parameter [hx = 0.9, 1.1, 1.3, from top to bottom). We can see 
that for 6ft- within la confidence interval from Maller et al. 2005, i.e. changing by 
over 40%, /3mle varies by no more than ~6% (for every given i^min)- What is more, 
for a fixed 6/^, there is a wide range of -R'mm for which /3mle is almost constant, as 
expected. For instance, if hx = 1.1, then 0.399 < /3mle < 0.408 for K„an e (4.6,6.1). 



85 



CHAPTER 5. OPTIMIZED ESTIMATION OF p PARAMETER 



rgff [Mpc/h] 

3.0 3.7 4.7 5.9 7.4 9.4 12 15 19 24 

0.52 I — I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I — I — I — I — I — I — I — I — I — I — I — I — I — 1 — I — 




0.38 I I 

3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 

Kmin [mag] 

Figure 5.3: Maximum-likelihood estimate of the /3 parameter for three values of = 
0.9, 1.1, 1.3, from top to bottom, as a function of increased minimum magnitude of the 2MASS 
sample, and for a constant maximum one, i^max = 13.6 mag. Note only slight variations of 
/3 for 4.6 < iTmin < 6.0. 

More generally, for parameter space bx G (0.9,1.3) and f^min ^ (4.6,6.1), we have 
0.393 < /3mle < 0.413. In order to account for possible systematic errors, such as 
those dealt with in Chapters 2 and 3, we take conservative Icr confidence intervals for 
our measurement and propose the following estimate of the (3 parameter based on this 
analysis: 

/3 = 0.40 ± 0.02 . [5341 

This value is slightly larger than, but consistent within the errors, with the one obtained 
from the analysis of the growth of the 2MASS dipole, presented in the previous Chapter. 
It is also in agreement with already mentioned results of the other authors who equally 
used 2MASS data: (3 = 0.40 ± 0.09 (Erdogdu et al., 2006a) and (3 = 0.35 ± 0.10 
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(NusSER ET AL., 2012), but has much better precision. We also note that repeating 
our analysis for '2MRS-like' data, i.e. setting the constant maximum magnitude of 
the sample at A'^ax = 11-75 mag, gives /3 = 0.43 ± 0.04 - a value consistent with 
those mentioned here, but with twice as big errors as our measurement from the whole 
sample. This increase in errors is caused by weaker correlation of the LG velocity and 
gravity if the maximum-magnitude cutoff is diminished to mimic the shallower extent 
of the 2MRS dataset, which results in less stable behavior of /3mle- 

If we now use our value of (3 to estimate the parameter = (/? we get 

from (5.34) 

= b]f^ X (0.19 ± 0.02) . {535] 

Finally, we can also constrain alone, using bx = 1.1 ±0.2 (Maller et al., 2005), 
to obtain 

= 0.225 ± 0.077 , {536] 

where, as in Chapter 4, such a big la error comes mostly from the uncertainty in the 
biasing parameter. This value of the density parameter is in good agreement with the 
estimates mentioned in Section 4.4. 

5.4 Summary and conclusions 

This Chapter presented a 'modern' approach to the measurement of the /3 parameter 
from the comparison of the LG velocity and acceleration. We used the same Two 
Micron All Sky Survey data as in Chapter 4 and took advantage of the fact that the 
growth of the 2MASS dipole is consistent with the ACDM model. This allowed us to 
focus only on the optimization of our measurement, for which we used the maximum- 
likelihood (MLE) approach. Application of this method for this problem is not a new 
idea; it has been used by various authors and in different variants for two decades now 
(Strauss et al., 1992; Schmoldt et al., 1999; Cieciel4g & Chodorowski, 
2004). However, only recently a really precise measurement of the /3 parameter has 
become possible, thanks to new datasets, such as the 2MASS XSC. 
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In our search for the optimal window, we followed the considerations of Chodor- 
OWSKI ET AL. (2008) and reduced the galaxy sample by removing the sources with 
the largest observed fluxes in the K band. We confirmed their predictions that the 
misalignment angle 9 between the LG velocity and the 2MASS dipole should decrease 
to a minimum value for a certain cutoff magnitude i^min- From our analysis we found 
that ^min = 9° ± 2° for -K'min = 6.3 mag. However, we have not managed to select 
a single optimal window for the MLE of /3. In particular, we have shown that the 
minimum of the misalignment angle is not observed for the same magnitude cut as 
is the maximum of the effective cross-correlation coefficient between the two vectors 
- contrarily to the theoretical expectations of COS, which were however strictly valid 
only in the absence of observational errors. The latter (estimated here via bootstrap 
resampling) are not very large, but still influence (slightly reduce) the cross-correlation 
coefficient. 

Although we do not select a single value of i^min at which the sample should be 
limited to obtain an optimal measurement, we can still determine the (5 parameter with 
excellent precision. Basing on the formula first derived by CC04 for the likelihood func- 
tion of cosmological parameters given v^^ and g-^^, we presented the equation for the 
maximum-likelihood estimate of /3 and argued that it depends very weakly on the bias- 
ing parameter h. By changing the observational window, we have observed very stable 
behavior of /3mle for a wide range of -ft'mm and hx- Our measurement is /3 = 0.40±0.02, 
which is consistent with the result presented in the previous Chapter and twice as pre- 
cise, and leads to the constraint on the density parameter Vt^^ = 0.23 ± 0.08. 

Such small errors in the /3 parameter measured via the ma:ximum-likelihood method 
are a great advancement when compared to older estimates, that used the data ob- 
tained by the IRAS satellite. Despite a rigorous statistical approach, neither Strauss 
ET AL. (1992), SCHMOLDT ET AL. (1999) nor even CC04 (who made the most pre- 
cise of these measurements) were able to reduce relative errorbars in /3 to less then 
~ 30%. The main reason was much sparser sampling of the IRAS dataset when com- 
pared to 2MASS XSC, as well as the fact that the galaxies shining bright in the far 
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infrared (wavelengths of IRAS) are much more biased tracers of stellar and total mass 
than those from the catalog we used. In fact, the choice of the near- infrared bands for 
the survey that later became 2MASS had been partly motivated by the fact that the 
distribution of its galaxies was expected to follow the one of the underlying density 
field, as apparently is the case (Westover, 2007; Davis et al., 2011). Adding to 
that the relatively dense sampling and unprecedented depth as for an all-sky survey, 
we have a catalog that is a great tool for precision cosmology, as we have shown in this 
and the preceding Chapter. 
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Conclusions and future prospects 



In this Thesis, we presented two approaches to measure the cosmological parameter 
P = Q'^^ /b, where Q^i is the mean matter density of the Universe and b is hnear 
biasing. Both our methods rehed on the comparison between the clustering dipole of 
the galaxies from the 2MASS Extended Source Catalog and the peculiar velocity of 
the Local Group. Before describing the two approaches and the results, we started in 
Chapter 2 by studying a particular non-linear effect, related to the fact that the nearby 
Local Void is partially hidden behind the Zone of Avoidance. We showed that even 
random filling of this region and elongation of the void along the line of sight should 
not bring crucial systematics into the determination of the amplitude and direction 
of the dipole. Next, in Chapter 3 we presented the 2MASS data that we used and 
explained in detail how we handled them for the purpose of the analyses that followed. 
Of particular importance was the statistical criterion to remove Milky Way sources 
from the sample, which - if kept - could seriously bias the measurement. 

Next two Chapters dealt with the application of the data. In Chapter 4 we stud- 
ied how the 2MASS clustering dipole grows when the sample depth is increased. By 
comparing the observations to the theoretical expectations of the ACDM cosmological 
model, we showed that once the proper window of the survey is included, the observed 
growth is consistent with the theoretical, conditional one (using the fact that the ve- 
locity of the Local Group is known independently). This comparison requires rescaling 
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the datapoints by the /3 parameter, which provided us with a straightforward method 
to estimate it. The result was /3 = 0.38 ± 0.04, where these la errors included such 
systematics as those dealt with in Chapters 2 and 3. 

Finally, in Chapter 5 we refined our measurement of /3. We used the maximum- 
likelihood machinery developed for this purpose first in the context of IRAS data by 
Strauss et al. (1992); Schmoldt et al. (1999) and Cieciel4G & Chodor- 
OWSKI (2004) and then adjusted for the 2MASS case by Chodorowski et al. 
(2008). We showed that removing from the sample the ~ 50 galaxies brighter than 
Kjnin = 6.3 mag allows to reduce the misalignment angle between the LG gravity 
and velocity to ~ 9°. However, the cross-correlation coefficient between the two vec- 
tors, even when corrected for observational errors, has its maximum for a significantly 
smaller J'^min, which precluded us from determining a single optimal window for the 
maximum-likelihood estimate. Still, for a wide range of minimum magnitude cutoffs 
we noted very stable behavior of /3 calculated within this framework, which allowed 
us to evaluate this parameter with double the precision of the method presented in 
Chapter 4. The final result is /3 = 0.40 ± 0.02, which has much smaller errorbars than 
not only the older maximum-likelihood measurements that used the IRAS data, but 
also when compared to more recent estimates for which 2MASS data were used, with 
various methods. 

Statistical precision (small errorbars) of a measurement is not necessarily tanta- 
mount to its accuracy, i.e. the ability to reproduce the real value of the estimated 
quantity. Unknown systematic errors can bias the estimate without the researcher(s) 
being aware of that fact. It seems however that our results are both precise and accu- 
rate, at least to a certain extent, as they are consistent with separate and independent 
analyses of a subset of the 2MASS catalog, namely the 2MASS Redshift Survey. Both 
Erdogdu et al. (2006a) and Nusser et al. (2012) found values of (3 that agree 
within the errors with those presented here. Additionally, our results are roughly com- 
patible with the findings of PiKE & HUDSON (2005) and Davis et al. (2011), where 
apart from the 2MASS data, also peculiar velocity catalogs where used to compare the 
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density and velocity fields. 

There is of course still room for improvement. First of all, it might be worthwhile 
to confront our findings with computer simulations. Generating thousands of mock 
'2MASS-like' catalogs and repeating the same measurements on them as we did here 
could help to validate our estimates of the errorbars, especially that simulated data 
would not suffer from such biases as the Zone of Avoidance or noise induced by Galactic 
sources that were not removed by our filter. Simulations could also be useful to verify 
how much our results may be biased due to the presence of nearby non-linear structures 
hidden behind the ZoA, other than the Local Void (such as the Norma Cluster, part 
of the Great Attractor, e.g. Woudt et al. 2008). 

Progress is also possible once new data become available. Firstly, an all-sky redshift 
catalog with a depth and coverage comparable to those of the 2MASS XSC would allow 
to mitigate the effects related to flux- limiting the dataset, by using the selection func- 
tion which one cannot calculate not having galaxy redshifts in hand. Unfortunately, 
there is currently no observational project aiming at creating such a catalog in the 
near future, i.e. within next several years. There may be however a way to overcome 
this problem, by estimating so-called photometric redshifts for all those galaxies in the 
2MASS XSC that do not have spectroscopic ones measured. The method bases on 
correlating redshifts with other properties of galaxies, such as colors in different bands, 
morphologies etc. The literature dealing with relevant methods is vast and is beyond 
the scope of this work. Anyway, a future photometric-redshift all-sky catalog based 
on 2MASS would be a great tool not only for such analyses as those presented here, 
but also for studies of galaxy bulk flows or the cosmography of the 'local' Universe, 
which still holds some secrets despite immense progress in extragalactic science we are 
now witnessing. We hope that more and more of these secrets will be revealed in the 
coming years. 
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A 

Data of the brightest 2MASS galaxies 



The following Table A.l presents the data of the 50 brightest galaxies in our sam- 
ple, beyond the Local Group. The list includes 2MASS XSC galaxies with Galactic 
latitude \b\ > 5° and the 3 brightest LGA objects, which are located in the ZoA. The 
parameters given are the following: (1) name of the galaxy; (2) Galactic longitude / and 
latitude b; (3) radial velocity in the Local Group frame and (4) redshift-independent 
distance estimate, both taken from the NASA/IPAC Extragalactic Database (NED); 
(5) morphological type, as given by SIMBAD; (6) extinction-corrected 20mag/sq." 
isophotal fiducial elliptical aperture magnitudes (J, H and Kg bands), without the 
offset of A = —0.2 mag that should be applied to compensate for lost flux (see Section 
3.2). Note that the distance estimate given here in column (4) is the 'mean' taken from 
NED, which is an unweighted average of several independent measurements that may 
vary significantly depending on the method. We thus do not quote the errors, which 
formally may be equal to several dozen per cent. For details see relevant NED entries. 
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APPENDIX A. DATA OF THE BRIGHTEST 2MASS GALAXIES 



Table A.l: First fifty 2MASS galaxies with the smallest Kg- 
band apparent magnitudes. Galaxies marked with an asterisk 
were added from the Large Galaxy Atlas. 



(1) 


(2) 




(3) 


(4) 


(5) 




(6) 




name 


in 




''^LGf [km/s] 


d [Mpc] 


type 


J [mag] 


H [mag] 


K [mag] 


NGC 253 


97.367 


-87.965 


276±3 


3.2 


Sc 


4.857 


4.132 


3.815 


M 81 


142.092 


40.9 


108±9 


3.7 


Sb 


4.782 


4.132 


3.897 


Centaurus A 


309.517 


19.418 


301±16 


3.7 


E 


4.927 


4.246 


3.947 


Maffei 1* 


135.862 


-0.551 


298±15 


3.5 


E 


4.791 


4.235 


4.128 


NGC 4945 


305.272 


13.34 


299±16 


4.0 


Scr 


5.503 


4.792 


4.468 


Maffei 2* 


136.498 


-0.326 


212±15 


3.2 


Sbc 


5.262 


4.721 


4.606 


M 82 


141.409 


40.567 


348±10 


3.9 


I 


5.738 


5.004 


4.633 


IC 342 


138.173 


10.58 


245±13 


3.1 


Scd 


5.479 


4.96 


4.674 


M 83 


314.584 


31.973 


301±13 


4.7 


SBc 


5.593 


4.952 


4.721 


Circinus* 


311.326 


-3.808 


189±15 


4.2 


Sb 


5.193 


4.701 


4.790 


M 104 


298.461 


51.149 


828±13 


10.5 


Sa 


5.896 


5.229 


4.99 


M 94 


123.362 


76.008 


352±3 


5.0 


Sb 


6.068 


5.408 


5.163 


M 64 


315.683 


84.423 


364±5 


5.4 


Sb 


6.299 


5.623 


5.381 


NGC 6946 


95.719 


11.673 


344±18 


6.5 


Scd 


6.138 


5.926 


5.42 


M 49 


286.923 


70.196 


873±10 


16.1 


E 


6.369 


5.732 


5.498 


M 51 


104.851 


68.561 


692±6 


8.0 


Sc 


6.486 


5.797 


5.588 


M 106 


138.319 


68.842 


507±5 


7.5 


Sbc 


6.498 


5.832 


5.592 


NGC 1316 


240.163 


-56.69 


1633±1 


19.0 


Sa 


6.546 


5.947 


5.68 


M 63 


105.997 


74.288 


546±4 


9.0 


Sbc 


6.681 


5.947 


5.722 


M 77 


172.104 


-51.934 


1191±4 


12.3 


S 


6.984 


6.267 


5.8 


M 60 


295.874 


74.318 


1017±8 


16.5 


SO 


6.739 


6.065 


5.815 


NGC 1269 


247.524 


-57.042 


703±8 


8.6 


SBa 


6.73 


6.094 


5.847 


NGC 3521 


255.532 


52.829 


602±12 


12.2 


Sbc 


6.779 


6.095 


5.847 


M 87 


283.778 


74.491 


1203±9 


16.7 


E 


6.806 


6.145 


5.896 



.continued on next page 
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(1) 


(2) 




(3) 


(4) 


(5) 




(6) 




name 


in 


Ml 




d [Mpc] 


type 


J [mag] 


H [mag] 


K [mag] 


NGC 3115 


247.783 


36.781 


422±15 


10.1 


SO 


6.808 


6.148 


5.92 


M 66 


241.962 


64.419 


590±9 


10.1 


Sb 


6.881 


6.196 


5.939 


M 101 


102.037 


59.771 


379±9 


7.2 


Sc 


6.893 


6.204 


5.94 


NGC 891 


140.382 


-17.415 


738±13 


10.2 


Sb 


7.277 


6.397 


5.97 


IC 356 


138.463 


13.11 


1105±13 


18.9 


Sbc 


6.891 


6.196 


5.99 


NGC 7331 


93.722 


-20.724 


1119±18 


14.2 


Sb 


7.075 


6.334 


6.08 


NGC 2903 


208.712 


44.54 


443±7 


9.4 


Sb/Sc 


7.002 


6.376 


6.089 


M 65 


241.332 


64.222 


671±9 


12.6 


Sa 


7.018 


6.343 


6.104 


NGC 4565 


230.761 


86.438 


1194±6 


13.3 


Sb 


7.224 


6.408 


6.124 


NGC 2841 


166.942 


44.151 


685±4 


17.8 


Sb 


7.083 


6.381 


6.152 


NGC 6744 


332.224 


-26.147 


706±8 


9.5 


S 


6.97 


6.434 


6.181 


NGC 3628 


240.852 


64.781 


709±8 


12.2 


Sbc 


7.351 


6.528 


6.211 


M 85 


267.713 


79.237 


649±5 


17.4 


SO 


7.144 


6.496 


6.249 


NGC 2974 


239.512 


35.014 


1691±19 


25.3 


E 


7.208 


6.847 


6.249 


M 86 


279.084 


74.637 


-348±8 


16.3 


SO 


7.164 


6.489 


6.272 


NGC 1023 


145.023 


-19.089 


828±12 


11.8 


SBO 


7.206 


6.547 


6.31 


NGC 2784 


251.971 


16.354 


402±39 


8.5 


SO 


7.205 


6.539 


6.317 


M 88 


282.33 


76.508 


2187±6 


19.5 


Sbc 


7.256 


6.604 


6.32 


M 84 


278.205 


74.478 


955±9 


17.4 


E 


7.224 


6.564 


6.332 


NGC 1553 


265.631 


-43.691 


869±17 


15.6 


SO 


7.245 


6.57 


6.345 


M 105 


233.49 


57.633 


762±9 


10.5 


E 


7.246 


6.573 


6.353 


NGC 2683 


190.455 


38.761 


365±3 


10.2 


Sb 


7.349 


6.645 


6.375 


NGC 4725 


295.085 


88.357 


1176±4 


13.6 


SBb/Sb 


7.341 


6.55 


6.38 


NGC 5195 


104.881 


68.488 


558±11 


(=M51) 


SBO 


7.329 


6.62 


6.389 


M 96 


234.435 


57.011 


744±10 


10.9 


Sab 


7.323 


6.634 


6.394 


NGC 1097 


226.915 


-64.68 


1198±5 


16.8 


SBbc 


7.372 


6.702 


6.429 



B 

Effective distance for a given flux 



In this Appendix we show how we calculate the effective distance of galaxies with a 
given flux S, knowing their luminosity function (LF) in the given band, ^{L), but not 
knowing their redshifts. Let us start by deriving the effective depth of a flux-limited 
sample. 

The number of galaxies dM in a volume element dV = 4 vr dr of a spherical 
sample limited by minimum flux 6*111111 (equivalent in the K band to some limiting 
magnitude -ft'max) is given by 

+00 

dN' = Anr^dr J ^{L)dL , (BA) 

where Lm\n = Anr minimum luminosity of galaxies in the sample at a 

given distance r. The mean depth of the sample is defined as 



/ r^dr 

which, for the Schechter LF (Schechter, 1976) with a faint-end slope a and char- 
acteristic luminosity L*, gives 

+00 

/ r3r(l + «,4 7rr2^iiiiii/L,)dr 

(R) = > 

/ r2r(l + a,4 7rr2 5iiiiii/i:,)dr 
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APPENDIX B. EFFECTIVE DISTANCE FOR A GIVEN FLUX 



where T{a, x) is the upper incomplete Gamma function. 

Using the distribution (B.l), we can also easily derive the median depth of a flux- 
limited sample, which we shall denote as R. We obtain it by solving the integral 
equation: 



dr = / — — dr 



dr 



dr 



R 



The following Table B.l presents the mean and median depth of the 2MRS (-ft'max = 
11.75 mag) and 2MASS XSC (-R'max = 13.6 mag) for the LF as given by Jones et al. 
(2006). The difference between the mean and the median is in that case of the order 
of a few percent. 



survey 


mean depth (R) 


median depth R 


2MRS 


104 Mpc/h 


96 Mpc/h 


2MASS XSC 


243 Mpc/h 


225 Mpc/h 



Table B.l: Mean and median effective depth of flux-limited surveys: the 2MASS Redshift 
Survey (2MRS, Kmax = 11-75 mag) and the 2MASS Extended Source Catalog (XSC, Kmax = 
13.6 mag), calculated from Eqs. (B.3) and (B.4), respectively, for the X-band luminosity 
function in the Schechter form with a = —1.16 and = —23.83 -|- 51og/i (Jones et al., 
2006). 



The effective depth of a flux-limited survey is however not a good measure of the 
effective distance of a thin shell in 'flux space' (in which the third coordinate is the 
flux, by analogy with the redshift space) that we seek. What we need is the mean 
value of distances of all galaxies with given flux S. The mean that we calculate will 
be thus a conditional one. This derivation is qualitatively the same as for the mean 
redshift of galaxies with a given flux, presented on pages 120-121 of Peebles (1993). 

We start by deriving the joint probability distribution of galaxy distances r and 
fluxes S. It is easily obtained by differentiating Eq. (B.l) with respect to luminosity 
L. Hence, the differential number S^N of galaxies with a LF $(-L) in a volume element 
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SV is given by 

5^N = 5LSV 

with SV = 47rr^ 6r. Now, passing from luminosity to flux, L = Anr'^S, for fixed r we 
have SL = Airr"^ dS. This gives the joint probability as 

fj2j\T 

p(r, S) = ^f— = 16 71^ $(47rr25). 

or Ob 

The conditional probability for r given S is 

vir\S) = ^ = ^6 7r^r^$(Wg) ^ 
' ^ p{S) /o+°"16 7r2r4$(47rr2^)dr ' 

hence the conditional mean for r given S will be 

+00 

/ $(47rr2S') dr 

(-)^ = is • 

J $(47rr2S') dr 


For the Schechter form of the luminosity function, let us define a characteristic dis- 
tance rs, which is the distance to a galaxy with a given flux S and the characteristic 
luminosity of the LF: 

Then for the Schechter LF, Eq. (B.8) simplifies to 

r(a + 3) ^ , 

/^\ o = r<^ — ^ — (B 10 J 

^r(a + 5/2) 

with r(a) being the Gamma function. 

Similarly, the conditional median for r given S, denoted here as fs, will be obtained 
from the implicit equation 

rs +00 

J ^{Anr'^S) dr = j ^{Aur'^S) dr , (bH) 

rs 

which for the Schechter LF is equivalent to solving 
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APPENDIX B. EFFECTIVE DISTANCE FOR A GIVEN FLUX 



with respect to f^, where 7(0, x) and T{a,x) are respectively the lower and upper 
incomplete Gamma functions. 

In the particular case of the K band, taking a = —1.16 (Jones et al., 2006), 
we obtain 



(r)5 = 1.056 rcj (bJs ) 

and 

fs = lM3rs. (baa ) 

Note that especially the median conditional distance gives a value very close to the 
characteristic distance rs, which could be the 'first-guess' effective distance (Tully 2008, 
private communication). 

For the purpose of the present work, we have decided to use the conditional median 
as our measure of the effective distance of galaxies with a given flux. This particular 
choice is somewhat arbitrary and influences the scaling of top axes in Figures 4.2, 4.4, 
4.6 and 4.7. However, since the difference between the conditional mean and median is 
very small (about 5%), this choice has a negligible impact on general conclusions and 
results of Chapter 4. 

Relating the flux to the K magnitude via Eq. (4.2), S = Sq 10~°'^^, and owing to 
the luminosity-magnitude relation = 47r(10 pc)^5'o 10~°'^^^*, we obtain 

reff = rs = 1-013 x io°-^(k-m,)-5 y^^^ ^ ( bIs) 

so finally, for = -23.83 + 5\ogh (Jones et al., 2006), 

rcff = 0.591 X 10°-^^ Mpc/h . (bJ6 ) 

For example, galaxies at the limit of the 2MASS catalog (with K = 13.6 mag) are 
assigned an effective distance of = 310 Mpc/h. 



102 



c 

Error estimation via bootstrap 

resampling 



Direct estimation of the errors of the measured dipole, as well as of the misalignment 
angle (Chapters 4 and 5), is very difficult if feasible at all. Many factors contribute to 
these errors, such as photometric and astrometric uncertainties in the catalog (which 
are however tiny), biases related to the data preparation procedure as described in 
Chapter 3, etc. For that reason we have chosen to use the widely adopted bootstrap 
resampling procedure to estimate the confidence intervals of our measurements. In 
general, bootstrapping is a statistical method for constraining the distribution of an 
estimator by sampling with replacement from the original sample, most often with 
the purpose of deriving robust estimates of standard errors and confidence intervals of 
population parameters (Efron & Tibshirani, 1994). In our case, it was the 2MASS 
catalog that was resampled. Based on the original catalog, we created 'mock' datasets 
by drawing with replacement from the original one as many sources as actually exist. 
Such a mock catalog thus created contained the same number of galaxies as the real 
2MASS XSC, however some of the original galaxies were missing, while some were 
repeated more than once, in line with the bootstrap ideology. For our purposes we 
created one thousand such catalogs (independent of each other), we checked however 
that the procedure converged already for 100 mocks. For each of these 1000 catalogs. 
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APPENDIX C. ERROR ESTIMATION VIA BOOTSTRAP RESAMPLING 



we applied the similar procedure as for the real data. Namely, for the purpose of 
the analysis presented in Chapter 4, we calculated the 2MASS clustering dipole and 
the misalignment angle by increasing the maximum magnitude of the sample. The 
values obtained for each of the mocks up to a given -R'max were averaged and used to 
calculate the la variance. These variances were applied as error estimates for the data 
in Figure 4.7 and in Subsection 4.3.3. Similarly, for the purpose of Chapter 5, we did 
the bootstrap for every subsample limited by a given -fTmin and -ft'max, which allowed 
us to estimate the errors in the measured gravity and (the cosine of) the misalignment 
angle (cf. Figures 5.1 and 5.2). 
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